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Abstract 



The electromagnetic and pionic transitions in mesons with heavy flavor charm (c) or bot- 
tom (b) quarks are calculated within the framework of the covariant Blankenbecler-Sugar 
(BSLT) equation. The magnetic dipole (Ml) transitions in the charmonium (cc) system 
are shown to be sensitive to the relativistic aspect of the spin-flip magnetic moment op- 
erator, and the Lorentz coupling structure of the QQ interaction. The observed rate for 
the Ml transition J/ip —^ rjd is shown to provide strong evidence for a scalar conflning 
interaction. On the other hand, the electric dipole (El) transitions are shown to be 
sensitive to the hyperfine splittings in the QQ system, and to require a nonperturbative 
treatment of the hyperflne components in the QQ interaction. 

In addition to the spin-flip Ml transitions, the single pion (tt) and dipion (tttt) widths 
are calculated for the heavy-light {Qq) D mesons, by employment of the pseudovector 
pion-quark coupling suggested by chiral perturbation theory. The pionic transitions 
D* — > DiT are shown to provide useful and constraining information on the pion-quark 
axial coupling g\. It is also shown that axial exchange charge contributions associated 
with the Qq interaction suppress the axial charge amplitude for pion emission by an order 
of magnitude. The models for tt and Ml transitions also make it possible to estimate the 
77-nucleon coupling from the transition D* — > DsTt'^, once the value of the tt^ — 77 mixing 
angle is known. 

Finally, the tttt dipion transition rates of the L = 1 D mesons are calculated, and are 
shown to make up a significant fraction of their total widths for strong decay. The tttt 
transitions between S'-wave charmonium (cc) and bottomonium (66) states are modeled 
in terms of a broad a meson or a glueball, with derivative couplings to pions. The effects 
of pion rescattering by the spectator quark are also investigated. 
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Short Introduction to the 
Papers 



I This paper presents a calculation of the Ml transition rates in the cc and bb systems 
within the framework of the nonrelativistic Schrodinger equation. A relativistic 
version of the single quark spin-flip magnetic moment operator is derived, along 
with the two-quark exchange current operators for Ml transitions. It is shown that 
the two-quark operator associated with a scalar confining interaction may provide, 
together with the relativistic single quark operator, a possible explanation of the 
empirically measured width of 1 keV for the transition J/tp ^ r]c'y. 

II The framework of the covariant Blankenbecler-Sugar (BSLT) equation is used to- 
gether with the pseudovector pion-quark coupling suggested by chiral perturbation 
theory to predict the widths for pionic transitions in the heavy-light (eg) D meson 
systems. It is found that useful and constraining information on the pion-quark 
axial coupling g"^ is provided by the D* Dir transitions. A satisfactory descrip- 
tion of the empirically measured pion widths of the L = 1 D2 meson is obtained. 
Also, the axial charge component of the amplitude for pion emission is shown to be 
suppressed by axial exchange charge contributions associated with the Qq interac- 
tion. 

III The chiral pseudovector Lagrangian, augmented by a Weinberg- Tomozawa term for 
dipion emission, is used to predict the widths for tttt transitions from the L = 1 D 
mesons. It is found that widths of several McV are expected for these transitions, 
in analogy with the experimentally well-studied decays of the strange K2 meson. 
It is thus expected that the tttt modes should constitute a significant fraction of the 
total widths of the L = 1 D mesons. 

IV The dipion transitions between 5-wave states in the charmonium (cc) and bot- 
tomonium [bb) systems are studied using a phenomenological model with deriva- 
tive couplings to pions. The dipions are modeled in terms of a broad a meson 
or a glueball. Effects of pion rcscattcring by the spectator quark are investigated 
and shown to be small for 2S IS transitions. The present experimental data 
on these transitions is shown to constrain the a meson mass to about 500 MeV. 
Finally, it is demonstrated that the anomalous double-peaked tttt spcctriim of the 
T(3S') T(IS') TTTT transition may be modeled in terms of a heavier ~ 1500 MeV 
scalar meson. 



VI 



V The empirically measured branching ratios for D* — > Dg 7r° and D* ^ Dgj are 
shown to provide a means of determining the strength of the rj coupling to quarks 
and baryons. This requires that the value of the tt" — 77 mixing angle is available, 
along with realistic models for the Ml and pionic transitions in heavy-light mesons. 
The value thus obtained for the 77-nucleon pseudovector coupling /^^r n is shown to 
be much smaller than that suggested by SU(3) symmetry, which is consistent with 
other recent phenomenological analyses. It is also shown that a significant Tj-charm 
coupling, if present, serves to increase the estimated value of f^NN- 

VI The two-quark exchange current operators that arise from the elimination of the 
negative energy components of the Bethe-Salpeter equation in the BSLT quasipo- 
tential reduction, are calculated for electromagnetic El and Ml transitions in heavy 
quarkonium systems. Although the exchange charge operators that contribute to 
El transitions are shown to be mostly negligible, the corresponding exchange cur- 
rent operators for Ml transitions are shown to be crucial, if agreement with the 
empirical width for J rjc 7 is to be achieved. This requires that the effective 
confining interaction couples as a Lorentz scalar, since an effective vector interac- 
tion is shown to yield a spin-flip magnetic moment operator only if the constituent 
quark masses are unequal. Consequently, in the system, the one-gluon exchange 
interaction also contributes a two-quark spin-flip operator. 
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Chapter 1 

Introduction 



It has been widely accepted, since the beginning of the 20th century, that the visible 
matter in the universe is composed of protons and neutrons (or baryons), and electrons 
(or leptons). However, the discovery of the positron in 1933, predicted by Dirac a few 
years earlier, suggested that short-lived, transient particles may exist alongside the stable 
protons and electrons. This was confirmed in 1936, when a heavier, unstable electron- 
like particle, the muon (/i), was discovered in cosmic ray experiments by Anderson and 
Neddermeyer. This discovery was followed up in 1947, when the existence of the neutral 
(7r°) and charged (tt^) pions, predicted earlier by Yukawa to be the carriers of the strong 
nuclear force, was confirmed by a similar experiment. These particles were the first 
ones of a large number of short-lived, unstable baryons, mesons and leptons which were 
subsequently produced in copious numbers by accelerator experiments. In particular, 
the pions were shown to be the lightest members of a new family of particles known as 
mesons, to denote that they are intermediate in mass between the baryons and leptons. 



1.1 The quark model 

Around 1960, the number of short-lived baryons (A, E ,A ,S ...) and mesons (tt ,K ,p ,ry ...) 
that had been discovered by accelerator experiments was overwhelming. This suggested 
that the hypothesis of Mendeleev could be extended to the baryons and mesons; Rather 
than being elementary, they might possess substructure and could perhaps be classified 
according to a "periodic table" of subatomic particles. This notion, originally put forward 
by Gell-Mann ^ and Ne'eman [2] among others, became known as the "Quark Model" 
and attempts to explain the observed properties (spin, isospin, electric charge, parity) of 
the mesons and baryons by arranging them into multiplets according to the symmetry 
group SU{3). It was found that three quark flavors, "up" (u), "down" (d) and "strange" 
(s) with spin 1/2 and fractional electric charges were required to accommodate all of the 
mesons and baryons known at that time. 

The experimental discovery 3 of the Q baryon, which was predicted by the quark model 
because of a gap in the "periodic table" of the baryons, soon provided dramatic confir- 
mation of the quark hypothesis. Although the quarks were at first only thought of as 
a useful theoretical tool, their actual existence inside the proton was confirmed by 
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deep inelastic e^p scattering (DIS) experiments at high energies. However, in spite of 
these remarkable successes, the quark model soon ran into a difficulty of symmetry. The 
spin-parity quantum numbers of the A resonance were empirically found to be consistent 
with a combined spin-flavor and configuration space wavefunction which is symmetric. 
This is inconsistent with Fermi-Dirac statistics, which requires that the total baryon 
wavefunction should be antisymmetric. 

This critical problem was finally circumvented by the introduction of a new property 
for the quarks, "color", which allows the wavefunction to be made antisymmetric by 
means of three color quantum numbers. In order to avoid an undesirable proliferation of 
unobserved states, a further constraint was placed, namely that the quarks only combine 
into colorless states (or singlet representations of the color SU (3)c group). This restricts 
the possible ways of combining quarks and antiquarks to hadrons, the simplest being 
qq (mesons), qqq (baryons) and qqq (antibaryons). Together with the proposal that 
eight spin 1 gauge fields, "gluons", should be associated with the new symmetry group 
SU(S)c, these notions were eventually developed into the theory of strong interactions, 
called Quantum Chromodynamics (QCD) 0. 

It was also realized that a fourth quark is required in the theory of weak interactions 
to explain e.g. the observed rate for the decay The fourth quark was 

eventually discovered in the form of narrow resonances |7j in November 1974 at center- 
of-mass energies of 3.1 GeV and 3.7 GeV in e'^e~ annihilation and, independently, in 
proton-proton collisions. These resonances, named J/ip and ip' , were interpreted as 
mesonic bound states of the new "charm" quark and its antiquark, cc. The charm quark 
turned out to have a mass of ~ 1500 MeV, and is thus much more massive than the 
~ 5 MeV u, d quarks and the ~ 100 MeV s quark. Later, as higher collision energies 
became available, an unexpected "bottom" (b) quark with a mass of ^ 4800 MeV was 
similarly discovered in the form of bb or T mesons. This again raised the question of a 
possible partner for the b quark, and indeed an extraordinarily heavy " top" (t) quark was 
finally detected [H] in 1995, by the proton-antiproton collider experiments at Fermilab. 
The t quark turned out to have a mass of 175 GeV, which makes it the most massive 
elementary particle known, and it is too short-lived for mesonic tt bound states to form. 

1.2 Quantum Chromodynamics 

In the theory of Quantum Chromodynamics (QCD), the interactions between quarks are 
mediated by eight massless vector bosons called gluons. However, a number of compli- 
cations effectively prevent the properties of hadrons to be predicted from QCD; First of 
all, the theory is nonlinear due to gluon self-interactions, and it describes systems that 
interact strongly enough so that perturbative methods are inapplicable. Only at the 
very highest energy scales, where quarks become asymptotically free and the coupling 
between them small, can the predictions of perturbative QCD be compared with exper- 
imental results. At low energies, the quarks interact strongly, are confined into hadronic 
bound states and acquire effective masses. These constituent quark masses are for the 
light u, d quarks of the order ~ 400 MeV. 

At present, the only way to analyze QCD at a fundamental level is the method of "lattice 
QCD" simulations, where the properties of hadrons are probed by means of numerical 
Monte Carlo algorithms. Although much progress is being made in the development 
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of more efficient algoritlims and tlie inclusion of dynamical fermions (unqucnchcd lat- 
tice QCD) into the simulations, the applicability of such methods is still limited by the 
huge demands on computing power. In such a situation, it is natural to attempt to 
understand the properties of hadrons by means of effective theories and phenomenolog- 
ical, QCD-motivated models. The physical motivation of such an approach is that the 
fundamental degrees of freedom of QCD are quarks and gluons, whereas low-energy ex- 
periments observe hadrons, which at least at long range interact by Yukawa-type meson 
exchange. It is, therefore, a reasonable expectation that the low-energy properties of 
QCD can be described in terms of an effective theory. In the limit of vanishing quark 
masses, QCD exhibits an invariance under chiral transformations that involve left- and 
right-handed quark fields separately. This symmetry is only approximate for quarks with 
a nonzero mass. The absence of parity doublets in the low-energy region of the hadron 
spectra suggests that this chiral symmetry is spontaneously broken at low energies |H] ■ 



1.3 Heavy flavor mesons 

Mesons that contain either two heavy quarks (cc, 66, cb) or one heavy quark and one 
light (eg, cs, bq, 6s) are special, since their masses lie in a region which is intermediate 
between the high-energy perturbative regime of QCD and the low-energy regime where 
the dynamics are governed by chiral symmetry breaking. Thus these heavy flavor mesons 
are likely to share features that are encountered in these two limits of QCD. One task at 
hand is then to determine phenomenologically, or from lattice QCD TITJ, the functional 
form, strength and Lorentz structure of the Qq and QQ interaction. 

Although the nonrelativistic Schrodinger framework can be applied to QQ systems 
with some success, a realistic treatment of the Qq system has a priori to be relativistic, as 
the velocity of the confined light constitutent quark is close to that of light. The papers 
presented in this thesis employ the covariant Blankenbecler-Sugar (BSLT) reduction [T^ 
of the Bethe-Salpeter equation, which has the advantage of formal similarity to the 
Schrodinger framework. An alternate approach is provided by the Gross quasipotential 
reduction which has been shown 14 to yield comparable results for the spectra of 
QQ and Qq mesons. 

However, as the mass spectra of the QQ and Qq mesons are well described JHl by a large 
number of phenomenological and QCD-motivated models, the spectrum alone cannot 
discriminate between different assumptions for the QQ and Qq interaction. Fortunately, 
as will be shown in this thesis, the observed rates for 7 and tt transitions in heavy 
flavor mesons may provide useful and constraining information on the quark-antiquark 
interaction, the quarkonium wave functions, and in particular, on the Lorentz structure 
of the effective confining interaction. As the negative energy components of the Bethe- 
Salpeter equation are eliminated in the BSLT (or Schrodinger) quasipotential reduction, 
two-quark transition operators that depend explicitly on the Lorentz structure of the 
QQ interaction appear as a consequence '16'. In particular, it will be demonstrated in 
this thesis that a pure scalar confining interaction compares favorably with the current 
empirical knowledge of Ml transitions in the charmonium system. It is noteworthy, 
that similar results have been obtained within the instantaneous approximation to the 
Bethe-Salpeter equation ^17 \ , which treats the negative energy components explicitly, i.e. 
without two-quark currents. 
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1.4 Transitions in heavy flavor mesons 

The transitions considered in this thesis include the radiative El and Ml transitions in 
the QQ systems, the Ml transitions in the heavy- light charm {D) and strange charm {Ds) 
mesons, and the tt and tttt transitions in the QQ and Qq systems. It is shown in papers I 
and VI that a possible solution to the long-standing overprediction [TH' by a factor ~ 3 
of the width for the Ml transition J /tp rjc'j emerges, if the two-quark exchange current 
operator associated with a scalar confining interaction is included along with a rclativistic 
treatment of the single quark spin- flip operator. 

On the other hand, the exchange charge contributions to the El transition rates 
are shown in paper VI to be highly suppressed by the large masses of the charm and 
bottom constituent quarks [201 ■ Similarly, the nonrelativistic predictions for the spin-flip 
Ml widths of the Qq mesons are shown in paper V to be unrealistic, as the confined 
light constituent quark requires a relativistic treatment. It is also shown that accidental 
cancellations in the single quark spin-flip operators render the Ml widths very sensitive 
to two-quark exchange current contributions. However, as the form of the Qq interaction 
is uncertain, the results are suggestive rather than definite, quantitative predictions. 

In the heavy-light D mesons, the excited states decay to the ground state predominantly 
through pion emission. In this thesis, the pionic transitions in the D mesons are described 
in terms of the chiral pseudovector Lagrangian which couples the pions to the light 
constituent quarks. It is shown in paper II that the D* — > Dt: transitions can provide 
useful and constraining information on the pion-quark axial coupling g^^. Also, the axial 
charge component of the amplitude for pion emission is shown to be highly affected by 
two-quark axial exchange charge contributions associated with the Qq interaction. The 
pionic transitions which are driven by the axial charge operator may, therefore, provide 
information on the Lorentz structure of the Qq interaction. In particular, it is shown that 
a scalar confining interaction has the effect of reducing the widths for such transitions. 

The chiral Lagrangian may, when augmented with a Weinberg- Tomozawa term for dipion 
emission, describe the tttt widths of the excited L — 1 D mesons. In this thesis the tttt 
widths of the D mesons are shown to be of significant magnitude compared to the widths 
for single pion emission. This is known to be the case for the strange K2 meson, where 
the empirical tttt width is ~ 1/2 of the vr width. This model for pseudoscalar emission 
has also been applied to the D* Dgir^ transition, which can then be used to extract 
the coupling of 77 mesons to quarks and baryons from the empirical branching ratios for 
those transitions, once an estimate for the tt^ — rj mixing angle is available. The value 
for the ?7-nucleon pseudovector coupling constant friNN so obtained, is shown to be much 
smaller than that suggested by naive SU{Z) symmetry, but consistent with other recent 
phenomenological analyses of e.g. photoproduction of 77 mesons on the nucleon. 

Whereas the dipion transitions in the D mesons may be modeled in terms of the chiral 
Lagrangian, the tttt transitions between 5'-wave cc or bb states are likely to involve a 
broad a meson or a glueball. It is shown, within a model where the coupling of dipions 
to heavy quarks is mediated by a broad and heavy scalar meson, that the empirical tttt 
energy spectra constrain the a meson mass to ^ 500 MeV. A possible explanation for the 
anomalous double-peaked tttt spectrum of the T(35') T(lS')7r7r transition is obtained, 
if the TTTT emission is described in terms of a heavier (~ 1500 MeV) scalar meson. 
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1.5 Notation and layout 

Throughout this thesis, the natural units with he = 1 and the (5^^ metric have been 
employed. The Euclidean (5^^ or Pauli metric assigns imaginary time components to 
four-vectors. The momentum four-vector k is thus of the form k = {k,iko), where the 
three-vector has been denoted by bold-faced type. However, for typesetting reasons, 
three-vectors in exponents have been denoted with an arrow, according to k. Also, in 
obvious cases the modulus |fe| has been denoted simply by A;. In the Pauli metric the 
square of a four- vector is of the form 



and the Dirac 7^ matrices are all hermitian with square equal to one. The explicit form 
of these matrices in the Pauli metric is then 7^ = (7,74) and 75 = 71727374, where 



Factors of i arc also introduced into the Dirac current and charge density operators 
to make them real- valued, and for Lagrangians which include a 75, in order to assure 
hermiticity. 

A number of abbreviations that are frequently used in this thesis are OGE (one-gluon ex- 
change), BSLT (Blankenbecler-Sugar-Logunov-Tavkhelidze), NRIA (non-relativistic im- 
pulse approximation) and RIA (relativistic impulse approximation). Excited states in 
the heavy quarkonium systems have been denoted either by the ^jj{nJ) or the primed 
notation, where the nth excited state is denoted by n primes, e.g. '0(35') = V"- Note 
that in the primed notation, the primes refer to radial excitations only. 

This thesis contains a summary which comprises six chapters, and reprints of selected 
research papers that have been signed by the author. Chapter 2 of the summary presents 
the Blankenbecler-Sugar quasipotential reduction, the QQ and Qq Hamiltonian models 
and the numerical results for the spectra of the heavy flavor mesons. Chapter 3 discusses 
the calculations of the electromagnetic El and Ml widths of papers I, V and VI, while 
chapter 4 presents the calculation of the pionic transitions in the D mesons of paper II 
and the estimation of the jj-nucleon pseudovector coupling fr^NN from paper V. Chapter 
5 deals with the tttt transitions in the D mesons (paper III) and the model for the tttt 
transitions in the QQ mesons from paper IV. Chapter 6 contains a concluding discussion. 



(1.1) 




(1.2) 



Chapter 2 



Models for the Spectra of QQ 
and Qq Mesons 

Although several phenomenological models that employ a nonrelativistic treatment of 
the heavy quarkonia have succeeded in describing many features of the cc and bb 
systems, a realistic treatment of the heavy-light mesons has a priori to be relativistic as 
the velocity of the confined light constituent quark is close to that of light. Also in the case 
of charmonium and bottomonium, the compact size of the QQ system causes the charm 
and bottom quarks to move with relativistic velocities, in spite of their large masses. 
The reason for this is the effective confining interaction, which has a string tension of 
~ 1 GeV/fm and confines the constituent quarks to a region of radius < 0.5 fm. In this 
situation, a quasipotential reduction of the relativistic Bethe-Salpeter equation suggests 
itself as a natural framework for a covariant description of the heavy quarkonium systems. 



2.1 The BSLT quasipotential reduction 

The field-theoretical scattering matrix S may be written in the form 

% = 5f,-i{2ii)H{Pf-P,)Mf^,, (2.1) 

where the second term on the r.h.s. has been defined, for notational convenience, with a 
minus sign. The scattering amplitude M is then defined as 

Mf^ = u{p'q)u{p'^) M u{pq)u{p^), (2.2) 

where pi and p'^ denote the initial and final four- momenta of the quarks. Note that the 
antiquark will be described throughout by positive energy spinors. The Bethe-Salpeter 
equation for the scattering amplitude M can then be written (schematically) in the form 

M = K + KGM, (2.3) 

or explicitly, for an arbitrary frame, as 

Mip',p,P)^K{p',p,P) + i f ^ Kip',k,P)G{k,P) M{k,p,P), (2.4) 



7 
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where P is the total four-momentum of the Qq system, and p, k and p' denote the initial, 
intermediate, and final relative four-momenta of the constituent quarks. In eq. 1)2. 4|l . K 
denotes the interaction kernel of the Bethe-Salpeter equation, which in the nonrelativistic 
limit corresponds to the potential defined for the Schrodinger equation. This can be seen 
by comparison of eq. H2.3|l and eq. (|2.1(l in the Born approximation. Also, G denotes the 
Green's function of the Bethe-Salpeter equation, which is here taken to be the free fermion 
propagator. When bound states are considered, the inhomogeneous term in eq. (|2.4|l 
is dropped. The second term of the Bethe-Salpeter scattering equation is illustrated, 
along with the choice of momentum variables for the Blankenbecler-Sugar quasipotential 
reduction, by Fig. 12.11 



W + A + k 









K 




m 



W - A-k 



W + p' 



W -p' 



Figure 2.1: Illustration of the choice of frame and variables for the derivation of the 
Blankenbecler-Sugar (BSLT) reduction of the Bethe-Salpeter scattering equation for 
unequal quark masses. The upper and lower quark lines are taken to have masses 
mq and mg, respectively. The four-momenta are defined as W = (0,iPo/2), A = 
(0, i[m'l - m^]/4:Wo), p = (p, ipo) and k = (fc, iko). 



It is instructive, in order to perform the BSLT quasipotential reduction, to introduce the 
variables presented in Fig. 12.11 and write the Bethe-Salpeter equation, schematically, as 
two coupled integral equations, 

M = U + UgM (2.5) 
U = K + K{G - g)U. (2.6) 

Here the quasipotential U is defined by eq. H2.6|l in terms of the Bethe-Salpeter propa- 
gator G and a three-dimensional propagator g. The propagator g is then constructed so 
that it has an identical elastic unitarity cut (right hand cut) as G in the physical region. 
The approximation U K will be employed here, in order to arrive at a major sim- 
plification of the Bethe-Salpeter problem. The propagators G and g will have identical 
discontinuities across the right hand cut if Disc G — 2i Im g. By means of the Cutkosky 
rules, Im g may then be obtained as 

27^2 

Img = ~——^[-f'^(W + k + A)+imQ][j''{W~k~A)+imq] 
(27r) 

<5(+) [{W + k + Af + m^] [{W-k-A)^ + m| ] , (2.7) 

where it has been indicated that only the positive energy roots of the arguments in the 
delta functions are to be included. The complete propagator g is then reconstructed by 
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means of the dispersion integral 

9 = - Im5(W^',A;,A'), (2.8) 

TT Jo q-p- le 

where W is defined as with go on sheU. Evaluation of the above integral yields the 
following form for the BSLT propagator g: 

1 J(fco) [l^EQjk) - 7O ■ fc - imq] [4Eq{k) + ■ k - ^m,-] 
^ 2 (27r)3 {EQ{k) + E^{k)) {k^ - - le) 

where the delta function ensures the condition fco = in the resulting three-dimensional 
integral equation. Note that the (in principle arbitrary) variable A was chosen so that this 
condition is realized also in the case of unequal constituent quark masses. By introduction 
of the positive energy projection operators A!^, the above propagator can be written in 
the form 

Sjko) 2mQmg Ag(fc)A^ (-fc) 
{2Trf EQ{k)+Eg{k) k^-p2 

This form is convenient when matrix elements are taken between positive-energy spinors 
according to 



X,V(p',p) = uq{p')u^{~p')M,U{p',p)uq{p)u^{-p), (2.11) 
which, together with cq. H2.5|l . yields the three-dimensional BSLT scattering equation 



(2J2) 

where V denotes the nonlocal interaction operator as obtained from the Feynman rules 
for Sfi using eq. (|2.1|l in the Born approximation. The above extension of the original 
BSLT equation to the case of unequal masses is similar to that of ref. 21 , which has 
been employed in ref. j22| for the case of AA^ scattering. 

The elimination of the negative energy components of the Bethe-Salpeter equation in 
the derivation of eq. (|2.12|l has been shown to give rise to two-quark exchange cur- 
rent operators that depend explicitly on the Lorentz structure of the quark- antiquark 
interaction. These may then contribute significantly to the strong and electromagnetic 
transition rates in the Qq and QQ systems. In particular, the exchange current oper- 
ator associated with the scalar confining interaction has been shown to be of decisive 
importance for the Ml transitions of heavy quarkonia (23]. It should be noted that the 
appearance of such two-quark operators depends on the type of quasipotential reduction. 
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Although eq. 12.1211 is a widely used quasipotential reduction of the type discussed in 
this thesis, there is in principle an infinite number of different ways to reduce the Bethe- 
Salpeter equation to a 3-dimensional form. Another commonly used reduction is the 
Thompson equation [23 j which differs from the BSLT equation by the choice of the 
dispersion integral H2.8|l . These have been shown to produce results that are very close 
to the full Bethe-Salpeter equation in ref. |2S]. There exists also a large variety of 
quasipotential reductions that differ in the choice of the propagator H2.10|l . which attempt 
to include the effects of intermediate negative energy states by various combinations of 
the negative energy projection operators |2f)j . 

It is also noteworthy that the Bethe-Salpeter equation in the ladder approximation has 
been shown [^T] , not to reduce to the desired one-body (Dirac) equation when one of the 
quarks becomes much heavier than the other. However, a large number of quasipotential 
reductions (e.g. Gross) are known to be closely related to the Dirac equation. This 
suggests that such reductions are more appropriate for two-quark systems with a large 
difference between the constituent masses, while the BSLT equation is ideal for quarkonia 
such as cc and bb. As the light constituent quarks in Qq mesons have masses that are 
lighter than those of the heavy quarks by factors of 3 — 10, then the Gross and BSLT 
reductions are expected to give results of similar quality, which indeed appears to be the 
case [Tl] . 



2.2 The BSLT and Lippmann-Schwinger equations 

As eq. H2.12|l is similar to the nonrelativistic Lippmann-Schwinger equation, except for 
the factor in parentheses, then it can be transformed into such an equation by means of 
the "minimal relativity" ansatz |2H1 

This yields the equation 

Tip',p)^V{p\p)- [ ^V{p',k) \ . T(fc,p), (2.15) 

which is formally identical to the Lippmann-Schwinger equation. Here /i stands for the 
usual reduced mass of the two-quark system. The advantage of eq. H2.15|l is that it can 
be transformed to a Schrodinger-type differential equation where the potential is given 
by eq. H2.14|l . This transformation gives the differential equation 

Ho - i>nim{r) = -V V'„im(r), (2.16) 

where -ffo is the kinetic energy operator of the nonrelativistic Schrodinger equation. The 
factor p2 can be expressed in terms of the total energy of the Qq state and the constituent 
quark masses mg and niq. 
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The eigenvalue e of the BSLT equation is obtained as 

where i? is the mass of the Qq state. The BSLT equation can thus be expressed as an 
eigenvalue equation of the form 

{Ho + Hint) ^nlm(r) = £ Vnlm(r-), (2.18) 

where the interaction Hamiltonian Hint is given in terms of the potential defined in 
eq. (I2.14|) . The introduction of the quadratic mass operator (|2.17|) leads to an effective 
weakening of the repulsive kinetic energy operator, which means that higher excited 
states will have lower masses in the BSLT equation than they would in the Schrodinger 
framework. The BSLT eigenvalue e, expressed in terms of the Schrodinger excitation 
energy ii'cx — E ~ (toi + m2), is of the form 



Ecx 

8^l 



(Ecx + 2(mQ + TO,)) (£^4 + 2£;cx(toq + "^g) + ^mqiriq) 



El^ + 2Ecx{mQ + rriq) + {mq + rnqY 



(2.19) 



where the expression in parentheses tends toward 8/i when toq, mq —f oo. This demon- 
strates that in the limit of heavy quark masses, or when the quark masses become large 
compared to the excitation energy E^x, the BSLT equation reduces to the nonrelativistic 
Schrodinger equation. 

Although the role of the BSLT potential V as given by eq. (|2.14|) is equivalent to that 
of the nonrelativistic, static potential in the Schrodinger framework, the multiplication 
of the full non-local interaction (in momentum space) by the minimal relativity square 
root factors is shown in the next section to have important consequences, not only for 
the numerical treatment of eq. H2.18|l but also for the modeling of the strong and electro- 
magnetic transitions between Qq and QQ states. In particular, the well-known problem 
of too singular and thus illegal hyperfine operators in the Schrodinger equation is shown 
to disappear in the BSLT framework. 



2.3 The QQ interaction in the BSLT framework 

The interaction between heavy quarks and heavy or light antiquarks is dominated by the 
(presumably linearly) rising confining interaction. The observed spectra of QQ mesons 
also require the presence of a short-range hyperfine interaction that gives rise to e.g. the 
J/'0— ?7c splitting. The one-gluon exchange (OGE) interaction [20] of perturbative QCD is 
a natural candidate for the QQ systems, whereas the origin of the hyperfine interaction in 
the Qq systems is less obvious. A recently suggested possibility is the pointlike instanton 
induced interaction proposed by ref. j30j . The interaction Hamiltonians used in this 
thesis in conjunction with the covariant BSLT equation are, therefore, of the form 

Hint = Konf + VoGE + Vinst , (2.20) 

with confining, OGE, and instanton induced components, respectively. The effective 
confining interaction is taken to have scalar Lorentz structure, while the OGE interaction 
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has vector coupling structure. In the nonrelativistic approximation, the effective linear 
confining interaction has the form (in the LS'-coupling scheme): 



cr 



mg + m? + rnqniq 



AmQrriqr 



r AmQUiq 



S-L+- "5 



{(Tq - (Tq) ■ L, 



(2.21) 



where the string tension c JSl is of the order ~ 1 GeV/fm. The above form contains 
also the momentum dependent terms from eq. H2.14|l up to second order in u^/c^. In the 
Schrodinger framework (i.e. without the minimal relativity factors), the numerical factor 
3/2 in front of the term would be 1. Likewise, the Darwin- Foldy term in eq. H2.21II 
would vanish. In addition to the familiar Thomas-precession term, an antisymmetric 
spin-orbit interaction also appears for unequal quark masses, which mixes the states 
with L = I and J = 1. 

The interaction components associated with the perturbative OGE interaction are, to 
order j (? in the nonrelativistic approximation, of the form 



Vqge — 




9 rnqraq 



where as denotes the strong coupling of perturbative QCD, and is the tensor operator 
'5i2 = 3((Tq • r){<Tq ■ r) — erg • cTq. In the Schrodinger framework, the coefficients for the 
contact and P^ terms would be — tt and 1, respectively. 

The instanton induced interaction, considered by ref. POI for systems with heavy quarks, 
consists of a spin-independent term as well as a ctq ■ CTq term which contributes to the 
pseudoscalar-vector splittings in heavy quarkonia. The effective instanton interaction 
derived in ref. [301 is of the form 



inst 



AA/g ^Mq 

An 



5{r) 



AM^"" ^Mq 

An 



S{r) (TQ ■ a-q 



9' 



(2.23) 



where the factors AMq and AMq denote the mass shifts of the heavy and light con- 
stituent quarks due to the instanton induced interaction. These shifts are, for light 
constituent quarks, of the order of the constituent quark mass (^ 400 MeV), and smaller 
(~ 100 MeV) for the charm quark. The parameter Mg"" controls the strength of the 
spin-spin interaction, which has the same sign as that from the perturbative OGE in- 
teraction. The parameter n represents the instanton density, which is typically assigned 
values around ^ 1 fm^^. The spin- independent term has scalar coupling for the light 
constituent quark line and a mixed scalar-70 vertex for the heavy quark. 
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2.4 Relativist ic QQ potentials 

For systems that contain light quarks, the above static interaction Hamiltonians have 
but qualitative value because of the slow convergence of the asymptotic expansion in 
v/c. Even for systems composed of heavy quarks only, the compact size of the wave 
functions lead to very large matrix elements in first order perturbation theory for the 
terms in eqs. H2.21|l and (|2.22|) . Therefore, it was chosen in ref. [201 to employ 
a local interaction model for the heavy quarkonium systems which takes into account 
the minimal relativity factors H2.14|l . as well as the relativistic effects due to the quark 
spinors and the running coupling of QCD. The central, spin-independent part of the 
OGE interaction is thus modified to 

VScEir) = -\- r dkj,{kr)WQ,'^^as{ei (2.24) 
where the following notation has been introduced for convenience: 

For the running QCD coupling as(fc^), the parameterization of ref . |31| : 

has been employed. Here the QCD scale parameter Aqcd and the dynamical gluon mass 
nig, which determines the low-momentum cutoff of the inverse logarithmic behavior of 
as have been determined by a fit to the experimental spectra of the QQ and Qq systems. 

In general, the relativistic effects in eq. (|2.24|l lead to a strong suppression of the short- 
range coulombic potential. On the other hand, the running coupling as, when employed 
according to eq. H2.26|l . increases the strength of the OGE interaction for large distances. 
The end result is, that the OGE interaction, when calculated using eqs. (|2.24|l and H2.26|l 
bears little or no resemblance to a coulombic potential, even for the heavy cc system. 
This can potentially have serious phenomenological consequences since models that em- 
ploy a short-range coulombic interaction have, in general, provided good descriptions 
of the cc and bb spectra. However, the spin-independent part of the instanton induced 
interaction H2.23|l has been shown in paper VI to provide the necessary short-range at- 
traction, even if the OGE interaction becomes weak. In principle, the effective confining 
interaction is also subject to similar relativistic effects, but in view of its long-range 
nature, their effect will be very small. 

The hyperfine components of the QQ interaction, as given by eqs. 1)2.21(1 and H2.22|l are 
usually treated as first order perturbations since their behavior for small r is too singular 
to allow for direct numerical treatment. Modification of those hyperfine components 
according to eq. (|2.14() leads to expressions, which are weaker and more well-behaved, 
and may consequently be fully taken into account. The employment of individual wave 
functions for each member in a given hyperfine multiplet was shown, in paper VI, to be 
important for a realistic description of several electromagnetic El and Ml transitions in 
the heavy quarkonium systems. 



(2.26) 
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The expressions for the local hyperfine components of the Qq interaction that take into 
account the minimal relativity factors and the running QCD coupling are, in configuration 
space, of the form 



TT r 7o k egeq 



eq + rriq cq + niQ 



(2.27) 
(2.28) 

^OGE = ^^Q- ^q r * 3o{kr) a,(fc2), (2.29) 

^OGE = ^Si2 r dkk^j2{kr)^as{k^), (2.30) 

for the spin-orbit, spin-spin and tensor components of the OGE interaction, and the spin- 
orbit (Thomas precession) term from the effective scalar confining interaction. Note that 
the expression l|2.28(l for the spin-orbit term associated with the linear scalar confining 
interaction is obtained by means of the representation — Snc/k'^ in momentum space. 
This can be understood as the Fourier transform of a modified linear potential cr e^^' in 
the limit A ^ 0. The integral ifT^ is convergent even if that limit is taken analytically. 
The above expressions are also free of singularities that require a perturbative treatment. 
If the QCD coupling as is taken to be constant, then the hyperfine components, as given 
by eqs. 12.27II - H2.30|) . reduce to the static expressions of eqs. (|2.21|l and 12.22|l for large 
distances. 

As the instanton induced interaction for Qq systems, as given by ref. j3()| . consists of delta 
functions, it has to be treated as a first order perturbation. Such a treatment is very 
unfortunate here since the repulsive kinetic energy as given by the BSLT quadratic mass 
operator H2.17|l is very sensitive to the ground state energy relative to the sum of the 
quark masses. A perturbative treatment of a strong attractive interaction component 
would thus effectively lead to unrealistically small level spacings between the higher 
excited states. In view of this, the delta function of eq. H2.23f) has been treated according 
to 

Mn. = -^^^^ rdkejoikr)WQ,^, (2.31) 

which effectively leads to a smeared-out form of the instanton induced interaction. While 
the presence of Wgq is naturally suggested by the BSLT minimal relativity factors, the 
niQ I eg factors are entirely phenomenological, and have been inserted to allow for better 
convergence of the above integral. In the limit of very large constituent masses (the 
static limit), the above equation reduces to the form (|2.23(l . The spin-spin component of 
the instanton induced interaction was found in ref. ,30, to be significant for the heavy- 
light Qq systems, but very weak for the heavy-heavy QQ mesons. Because of this, the 
spectra shown in Table [TTI and Fig. l2.2l do not include that interaction. The calculated 
Qq spectra employed in paper II, shown for the D meson in Fig. 12.31 do not include 
the instanton induced interaction, since sufficient attraction was provided there by the 
OGE interaction, although at the price of an unrealistically large value for the QCD scale 
parameter Aqcd- 
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2.5 Spectra of heavy flavor mesons 

The cc, bb and spectra that are shown in Table \TJ\ and Fig. 12.21 have been obtained 
by solution of the BSLT equation for a linear scalar confining interaction, and OGE + 
instanton components modeled according to the expressions given in this section. The 
hyperfine components have been taken fully into account, so that all the states given 
in Table 12.11 are represented by different radial wave functions. This model has been 
employed for the calculations of the electromagnetic transitions in paper VI as well as 
the dipion transitions in paper IV. 

Table 2.1: Calculated and experimental cc, bb and states rounded to the nearest 
MeV, as obtained in papers IV and VI. The states are classified according to excitation 
number n, total spin S', total orbital angular momentum L and total angular momentum 
J. The experimental values are from ref. |32j, except for the recently observed |^ ric{2S). 
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Exp (cc) 
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4015 
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9458 


9460 


3099 


3097 


6361 


2^ Si 


10030 


10023 


3678 
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3^ Si 


10377 


10355 


4040 


4040 ± 10 


7244 


4^ Si 
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10580 


4319 


4159 ±20 ? 


7500 
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3513 




6754 
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10266 




3912 




7126 


3iPi 


10552 




4211 




7401 
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9855 


9860 


3464 


3415 


6723 
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10244 


10232 


3884 




7107 


S^Po 


10535 




4192 




7387 


l^Pi 


9883 


9893 


3513 


3511 


6751 


2^ Pi 


10263 


10255 


3913 




7125 


3^ Pi 


10550 




4213 




7400 


l^Pa 


9903 


9913 


3540 


3556 


6770 


2 3P2 


10277 


10269 


3930 




7136 


3 3P2 


10561 




4226 




7410 




10158 




3790 




7009 




10149 




3784 




7006 


l^Di 


10139 




3768 


3770 ± 2.5 


6998 



Although preliminary, the measured mass of the B^ was reported in ref. [31| as 6.40±0.39 
GeV, which is about ^ 100 MeV higher than the predicted 6308 MeV, and most other 
models ;35 give even lower masses for the B^ ground state. However, the predicted B^ 
spectrum agrees very well with the QCD-inspired model of ref. . 
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The quality of the calculated QQ spectra in Tabic is generally quite satisfactory, as 
both the ip' — J/ip and J/ip — rjc splittings are given realistically. In particular, the •qd'^S) 
state has recently been reported by the BELLE collaboration j33j with a mass of about 
3650 MeV. This suggests that the spin-spin splitting is much smaller for the 25* states 
than for the J/ijj and the r^c, a feature which is well described by the present model. 
The main difficulty is the prediction of the hyperfine splittings in the L = 1 multiplet of 
charmonium. Table im and Fig. 12.21 indicate that the splittings are underpredicted for 
cc but in reasonable agreement with experiment for hh. This problem can be traced, in 
part, to the weakness of the OGE tensor interaction as given by eq. (|2.30(l . 



A 



Mb 
Me 

QCD 



An 

(AM,,: 



Paper VI 


Other models 


4885 McV 


4870 MeV [Tl| 


1500 MeV 


1530 MeV HH 


260 MeV 


200-300 MeV [S] 


290 MeV 


rug > Aqcd eh 


890 MeV/fm 


912 MeV/fm [m 


0.084 fm^ 


- 0.05 fm^ inni 


0.004 fm^ 


? 



Table 2.2: Quark masses and coupHng 
constants used for the calculated spec- 
tra in Fig. 12.21 The values should be 
considered as best fits within the BSLT 
model to the empirical cc and hh spec- 
tra. 



The heavy quark masses are close to those obtained by Roberts et al. in ref. within 
the framework of the Gross equation. The values of Aqcd and nig are in line with those 
suggested in ref. [3l] j while the string tension c is somewhat smaller than that suggested 
by the lattice QCD calculations of ref. ;37 • The strength of the instanton induced 
interaction in the cc system is comparable to the estimate given in ref. |3()| . In spite of 
the generally satisfactory results, perfect agreement with experiment had to be sacrificed 
in the hh system in order to obtain an optimal description of the cc spectrum with the 
same set of parameters. As this nevertheless is a small effect, the results indicate that a 
flavor-independent confining interaction is a reasonable first approximation, in contrast 
to the instanton induced interaction, the strength of which depends explicitly on the 
quark flavors involved. 

The spectra of the heavy-light Qq mesons that were obtained in ref. within the 
framework of the BSLT equation have been used here for calculation of the Ml and pion 
widths of the Qq states. The D meson spectrum so obtained is shown in Fig. 12.31 Al- 
though quite satisfactory agreement with the empirical Qq spectra was achieved, this was 
only at the price of a very strong OGE interaction and the introduction of a negative 
constant into the scalar confining interaction, which was treated as a free parameter. 
The hyperfine components of the OGE and scalar confining interactions were treated 
as perturbations, while the instanton induced interaction was dropped since the OGE 
interaction was found to give sufficient attraction to account for the empirical spectra. 
However, later (unpublished) calculations of the Qq spectra have indicated that a non- 
perturbative treatment of the OGE spin-spin interaction strongly favors the inclusion 
of the instanton induced spin-spin interaction proposed by ref. |30j . in which case the 
spin-independent term of that interaction may also lead to a more realistic strength and 
low-momentum behavior of a,. 
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Figure 2.2: Calculated and experimental bb and cc spectra. All states are given in 
MeV, and correspond to the data in Table ITTl The thickness of the lines denoting the 
experimentally determined states indicates the uncertainty in the mass of the state. Note 
that the identification of the 4^5*1 state in charmonium is uncertain, and may actually 
be a 2^Di state, or a mixture of the two. The threshold for DD decay is at ~ 3750 
MeV, and for BB decay at ~ 10500 MeV. 



The small number of empirically known Qq states makes a determination of the quality 
of a given model difficult. However, the most significant unsolved problem in the case 
of the Qq spectra is the ordering of the L = 1 multiplet, of which only two resonances 
have been detected so far. These probably correspond to spin-triplet states with J = 2 
and J = 1 {jq = 3/2 in the heavy quark limit), and are denoted D2 and Z?i, respectively 
(note the notational confusion) . The empirical fact that the D2 is higher in mass by ~ 40 
MeV then suggests that the ordering of the L = 1 states in the D mesons is similar to 
that observed for the cc and bb systems. It is reassuring that this result is consistent with 
lattice QCD calculations of the spin-orbit splittings in heavy-light mesons [221 although 
the problem of poor convergence of such calculations is still not solved to satisfaction. In 
view of these results, it appears that the possibility |40| of spin-orbit inversion in heavy- 
light mesons is not realized. It should be noted, however, that the spin-orbit splittings 
as calculated from static expressions like eqs. (I2.21|) and H2.22() are unrealistic because of 
the low masses of the light constituent quarks. An unphysical dominance of the Thomas 
precession associated with the scalar confining interaction may then suggest that the 
spin-orbit splittings of the L — I Qq states should be inverted. 
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Empirical and caicuialed D meson spectra (Slate energies in IvleV) 



Experimental ■ 
Caicuiated 



Figure 2.3: Experimental 
and calculated D meson 
spectra from ref. PH|- It 
should be noted that the ex- 
cited 25' state D** at ~ 2630 
MeV, which was reported by 
ref. was not detected by 
a subsequent search 021 and 
may therefore not exist at 
that energy. 



It is instructive to compare the parameters of ref. i38 with the values suggested by the 
more realistic calculation (paper VI) of the QQ spectra in Table IT^ The heavy quark 
masses nic — 1580 MeV and mi, — 4885 MeV of ref. though slightly higher, agree 
quite well with those of Table 12.21 On the other hand, the masses of the light quarks 
were obtained as TO„,d = 450 MeV and nis = 560 MeV, respectively. Although the light 
quark mass is higher than the usual phenomenological value of ~ 350 MeV, it is close to 
the value 420 MeV which has been employed [30] for the instanton induced interaction. 
It has also been suggested 031 that a natural value of the constituent quark mass is one 
third of the A mass, rather than the nucleon mass. 

The parameters Aqcd and nig were obtained as 280 MeV and 240 MeV in ref. [55] 
and lead to a much stronger coupling as than is necessary in paper VI. Most other 
analyses suggest a much weaker form. The phenomenological consequences of a 
very strong OGE interaction for light constituent quarks are potentially serious [HI, since 
incorrect ordering of the positive and negative parity nucleon and A states may result. 
Moreover, the relativistic damping of the short-range part of the attractive coulombic 
OGE potential in ref. is such that the OGE interaction alone cannot describe e.g. 
the IS — IP level splittings in the Qq spectra. It was therefore necessary to push the 
other parameters of that model to their limits. 

In this situation, the instanton induced interaction of ref. [30] suggests itself naturally as 
it provides both a strong attraction in the 5- wave and contributes to the D* —D splitting. 
Also, the OGE potential of eq. H2.24|l for a system composed of two light (~ 400 MeV) 
quarks, becomes depleted for distances equal to the meson radius. The central OGE 
component may therefore only play a minor role in the dynamics of light constituent 
quarks, while for the Qq mesons, the short-range attraction may turn out to be best 
described by a combination of OGE and instanton induced components. 



Chapter 3 



Electromagnetic Transitions 



The radiative transitions in heavy quarkonium (cc, bb, cb) systems have drawn much 
theoretical interest ISIEH], as they can provide direct information on both the heavy 
quarkonium wave functions and the QQ interaction. As reasonably reliable empirical data 
now exists for a number of transitions in both the cc and bb systems |46| . a fair assessment 
of the quality of theoretical models is already possible. The measured 7 transitions in 
the charmonium (cc) system include the El transitions Xcj J /''Pi and -0' — > XcJl^ as 
well as the spin-flip Ml transitions J /ip —> rjd and "0 ' rjd. The situation concerning 
the analogous transitions in the bottomonium (bb) system is, however, less satisfactory as 
the total widths of the Xbj states are not known, and none of the spin-flip Ml transitions 
observed. 

Previous calculations of the El widths of heavy quarkonia have demonstrated that the 
El approximation leads to overpredictions of most transition rates, and that this over- 
prediction can be significantly reduced, if not entirely eliminated, by the consideration of 
relativistic effects. On the other hand, theoretical predictions for the Ml transitions have 
remained unsatisfactory for a long time ^Hl as the width for J/ip ^ t]c1 has typically 
been overpredicted by a factor ~ 3. However, calculations of Ml widths using the non- 
relativistic Schrodinger equation in ref. |47j and paper I have demonstrated that the Ml 
transitions in charmonium are sensitive both to the relativistic aspects of the spin-flip 
operator as well as the Lorentz structure of the QQ interaction. The results obtained 
in papers I and VI suggest that a scalar confining interaction may explain the observed 
width of ~ 1 keV for J/ip ^ rjcj, provided that an unapproximated single quark spin-flip 
operator is used. This conclusion is supported by the calculation of ref. ^7], which is 
based on the instantaneous approximation to the Bethe-Salpeter equation, even though a 
quantitative understanding of the radiative transitions in charmonium was not achieved. 

The situation concerning the Ml transitions in the heavy-light Qq mesons is more un- 
certain because of the scarcity of reliable empirical data. Only recently has a first mea- 
surement of the width of the D* state been published, which allows a determination of 
the partial width for the Ml transition D^* D^j from its known branching fraction. 
Even though the total width of the is not known, the partial width for Z?°7 
can be inferred from the measured width of the D^* and model calculations of the pionic 
widths of the D* mesons, as the relative branching fractions for tt and 7 emission are 
known. 
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As the velocity of the hght constituent quark in the Qq systems is close to that of light, the 
relativistic corrections to both single quark and two-quark operators will a priori be large. 
It is shown in paper V that the large relativistic corrections to the single quark spin-flip 
operators, that yield unfavorable results for the Ml widths of Qq mesons, are in general 
counteracted by the two-quark operators associated with an interaction Hamiltonian that 
consists of OGE -I- scalar confinement components, thus allowing for better agreement 
with experiment. It is also suggested that the instanton induced interaction for Qq 
systems of ref. [SU] may have a favorable effect on the predictions for the Ml widths of 
the D* mesons. 



3.1 Charge density and electric dipole operators 

The electromagnetic transition amplitude for a two-quark system, in the impulse approx- 
imation, is of the form 

T/. = -1 Sr,d\2 ^}iri,r2) e- [e^«"^"^ Ji (9) + e'^'^^^ ^2(9)] ^.(r-i,r2), (3.1) 

where q and i denote the momentum and polarization of the emitted photon, respec- 
tively, while ipi and iff denote the orbital wave functions of the initial and final heavy 
quarkonium states. In the above equation, Ji and J2 denote the single quark current op- 
erators of quarks 1 and 2, respectively. By Fourier transformation, the current operators 
may be rewritten as 



Jiq) = / dVe'«^-^~ j(r') (3.2) 
d^r'r'ij ■ V) e'^-^' ~ [ d^r'e'^'^' r'{\/ • j), (3.3) 



from which the El approximation is obtained if the exponentials in eq. H3.3(l are dropped 
(i.e. q 0). Application of the continuity equation then gives V • j = iujp. For nonzero 
q, the second term in eq. H3.3|l has to be retained without approximation. Note that 
J (q) is taken to contain the quantity in square brackets in eq. (|3.1() . Application of 
eq. 13.3|l together with eq. (|3.1|l then leads to the following form for the amplitude of a 
7 transition, 

Tfi ^ i\q\ J d^rid^r2(p*f{ri,r2) i ■ d{ri,r2)<fii{ri,r2). (3.4) 

The dipole operator d(ri,r2), 

d(ri,r2)= y"dVV«-^^V'p(r',ri,r2), (3.5) 

reduces to the El approximation in the limit q — )■ 0. In general, the charge density 
operator p{r') contains, in addition to the single quark contribution psq, an exchange 
part Pox, which arises from the two-quark currents that are illustrated by the diagrams 
in Fig. 13.11 A necessary constraint is that two-quark contributions to the charge density 
must have vanishing volume integrals. The dipole operator that corresponds to the single 
quark charge operator psq{r' ,r) = pi{r',ri) + P2{i'',r2) may be expressed as 

^ d^r' e^'^f'" r' e^^'^^^^-^^') + p2{q) r' e^'T-i-^-'"^^ , (3.6) 
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Figure 3.1: Two-quark contributions to the QQ current and charge density operators. In 
the decomposition of the QQ^ vertex, the irreducible two-quark contributions give rise 
to exchange current operators, the most important of which are iUustrated by the upper 
Born diagrams. In order to obtain the correct two-quark contribution, the positive energy 
part of the intermediate propagator is subtracted in the lower Born diagrams, since that 
part is already accounted for by the impulse approximation. Note that similar diagrams 
describe photon emission by the heavy antiquark. In the case of the QQ interaction, 
the scalar confining and vector OGE components have been taken into account. The 
contributions from the instanton induced interaction have not been considered but are 
in any case small for the QQ systems. 



which upon evaluation yields 

dsq(ri, ra) = lim [n e'^'^' pi{q) + e'^^-^i iV,-pi(q)] + (1 ^ 2), (3.7) 

where qf refers to the physical photon momentum of each transition. The above form 
reduces to the El expression by the substitution qj 0. The nonrelativistic single quark 
dipole operator is therefore of the form 



dsqiri,r2) = Qi ri e 



+ Q2r2 e^«^ '-% 



(3.8) 



where Qi is the charge of the heavy quark, while Q2 denotes that of the heavy antiquark. 
Insertion of eq. H3.8|l into eq. H3.4|l yields the single quark dipole operator 



dsqir) 



Qim2 - Q2'mi 
mi + m2 



r e 



lqfr/2 



(3.9) 
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The charge density operator in eq. 13.7|l is, to second order in w/c, of the form 



Psq — Qi 



1 - 



8rn^ 



ZCTl ■ p\ X Pi 

Am? 



(1 



(3.10) 



where the second term on the r.h.s. is the relativistic Darwin-Foldy term. The effect of 
this term is very small because of the large masses of the heavy constituent quarks. The 
justification of this expansion hes in the small coefficient of that term; It has been shown 
e.g. in papers I and VI that such an expansion cannot be used for the magnetic moment 
operator. 

If the two-quark exchange charge operators from the Born diagrams in Fig. 13.11 are 
decomposed as Pcxir' ,ri,r2) = Pexiir' ,ri) + Pex2('"', ''2), then the contribution from 
quark 1 may be expressed as 



Pcxiir',ri) = 



(27r)3 



fk2_ 
(27r)3 



Pcxi(g,fe2). 



(3.11) 



Here fe2 is the momentum transfered to the heavy antiquark and r is defined as ri — r2. 
The exchange charge contribution to the two-quark dipole operator 



'icx(r-i,r2) - / dVe'?-/-*''rVcx(r',ri,r2) 



(3.12) 



from quark 1 may then be expressed as 

d3fc2 



dexi (ri) = ri e 



iqf-ri 



(27r)3 

■ lim 

1^1} 



e~^''^->exl(g/,fc2) 



(27r)3 



'/Ocxi(q', ^2) 



(3.13) 



which again reduces to the El approximation by the substitution qf — )■ 0. 



The exchange charge density operators that are associated with the QQ interaction have 
been extracted in ref. |19j . for different Lorentz invariants for equal-mass systems. When 
generalized to unequal quark masses, the required operators are obtained as 



Pox = 



4™^ 

Qi 

Ami 



q ■ k2 2 q ■ k2 (Ti ■ (72 

3 



mi 



TO2 



Vgik2 



Q2 

4m| 



q ■ ki 



m2 

2 g ■ fcl (Tl ■ <T2 

3 mi 



(3.14) 
(3.15) 



Vg{kl). 



In the above expressions, Vc and Vg denote the Fourier transforms of the confining and 
OGE interactions, respectively. Evaluation of eq. H3.13|l thus yields the dipole operators 



d 



Conf 



(r) 



m2 



mi 



Am^ mi + m2 



Q2 



AmS mi -f m2 



r VrXr) e'^r^l-^ , 



Qi 


C - 


2cri..2\ 


Q2 f 1 


2 CTi - (T2 


Ami 


\mi 


3 m2 / 


4m| \m2 


3 mi 



r e 



iqf-r jl 



(3.16) 

r dr 



Here Vc{r) is the linear confining interaction, while Vg{r) denotes the form of the OGE 
interaction in configuration space, which is taken to include the effects of the running 
couphng of QCD. 
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3.2 Current density and magnetic moment operators 



In the impulse approximation, the spin-flip magnetic moment operator for Ml transitions 
between S'-wave heavy quarkonium states have been derived from the amphtude 

Tf, = -i27rf S^Pf-P^^qf) f dVv,}(r)£. [e^^^^/^ ^^(q) + e^^^-Z^ ^^(q)] (3.17) 



where r = ri — r2. Expansion of the exponential in eq. (|3.2|l according to ~ 1 + • r' 
then yields the Ml and E2 amplitudes for photon emission. Upon isolation of the Ml 
contribution, the matrix clement for J/ip rjcj and ^p' — > rjc'y may be written in the 
form 

Mf^ = i J d^r ip}{r) q x e ■ fi^^ ipi{r), (3.18) 
where fi^f denotes the spin-flip part of the magnetic moment operator 



fJ'^\ J d^r' r' xj{r'). 



(3.19) 



In eq. H3.19(l . the current operator consists of a single quark contribution j and a 
two-quark contribution j^-^, which arises from the pair terms given in Fig. 13.11 The 
corresponding single quark magnetic moment operator may be expressed as 



which yields 



d\ 
(2^ 



fJ-s 



d\' 



lim 



7q X [t 



iq-r/2 



(q)) 



(3.20) 



(3.21) 



The magnetic moment operator is given by eq. (|3.21() in the nonrclativistic impulse 
approximation. However, previous work has demonstrated that the static magnetic mo- 
ment operators of the baryons receive large corrections from the canonical boosts of 
the constituent quark spinors |48| . Furthermore, it has been shown in paper I that the 
nonrclativistic impulse approximation does not provide a satisfactory description of the 
spin-flip magnetic moment operators for QQ systems, even though the masses of the 
charm and bottom constituent quarks are large. The matrix element that corresponds 
to eq. I|3.18(l in the relativistic impulse approximation is of the form 



d^P 
(2^ 



(3.22) 



where the final and initial state coordinates r' and r are defined as r[ — r'2 and r\ — 
respectively. In eq. H3.22(l . the momentum variable P is defined as P = (p'-|-p)/2, where 
p' and p are the relative momenta in the representation = Pil2 + p,p2=Pi/1 — p 
and p'l = P f /2 + p\ p'2 — P j /2 — p' . The relativistic single quark magnetic moment 
operator that appears in the matrix element H3.22|l is of the form 



-2V,x 



(ji(q,P)+J2(q,P)) 



(3.23) 



where the single quark current operators (q, P) are now treated without approximation. 
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In the nonrelativistic approximation, the spin-dependent part of the single quark current 
operator is given by 



f sq 



— (cTi +0-2) X q 



Qi , Q2 



2toi 2m2 



— (cTi - (T2) X q 



Qi Q2 



2mi 2m2 



(3.24) 



where the first term describes the magnetic moment of the two-quark system whereas the 
second term is the spin-flip operator for an Ml transition in the nonrelativistic impulse 
approximation (NRIA). In order to obtain the relativistic single quark current operator 
to be used with eq. (|3.23|l . the nonrelativistic current operator for quark 1 should be 
replaced according to 



3i 



i(Ti X {p[-Pi)] 




(3.25) 



Pi 



E[ + nil El + nil _ 

and similarly for quark 2. In the above equation, the energy factors of the quarks are 
defined as Ei = \/pi + and E'l = yPi^rml ■ The spin-flip magnetic moment oper- 
ator in the non-relativistic impulse approximation (NRIA) may be obtained by insertion 
of eq. (|3.24|l into eq. (|3.21|) . giving 



A* sq 



Qi 

2mi 



Q2 
2m2 



(cTi - (T2). 



(3.26) 



The corresponding operator in the relativistic impulse approximation (RIA) has been 
considered in refs. [^HlEEl) and may for transitions between 5- wave states be expressed 
as 



Rol 
A* sq 



2m 



Q2 f-y 
J2 



2m 



(cTi - 0-2), 



(3.27) 



where the relativistic factors /7 are defined as 



J.7 ""z 

•'^ " 3^ 



(3.28) 



where Ei denotes the energy factor Ei ~ y P + m.f . This shows that a relativistic 
treatment will lead to an effective weakening of the NRIA spin-flip operator. 

In addition to the above single quark current operators, the pair terms in Fig. 13.11 also 
give large contributions to the magnetic moment operators of mesons and baryons 
However, in the case of the magnetic moments of the baryons, additional complications 
are known to arise from flavor dependent meson exchange interactions which also con- 
tribute signiflcant exchange current operators 48 . If the exchange current operators of 
Fig. 13. H are decomposed as Jcx(9j ki,k2) — JcxiC^i ^2)+3c^2{Q^ ^1)^ then the contribution 
to the two-quark magnetic moment operator may be written in the form 



Me 



(27r)6 



^iq-{r/2-r') ^ik-i-r 



r'xj^^^iqM) + (1^2) 



(3.29) 
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where it is again understood that r ^ — r in the contribution from quark 2. Evaluation 
of the above equation leads to an expression analogous to eq. H3.21|l . 



— lim 



iq-r jl 



(27r)3 



, p-«9-r/2 / d^kl ^iki-r 



(3.30) 



As the exchange current operators for most Lorentz invariants do not depend explicitly on 
the photon momentum q, one notable exception being that for the scalar invariant |49| . 
then the exchange magnetic moment operators turn out to be difficult to calculate directly 
from eq. H3.3()|l . It has therefore been shown in paper VI that a convenient way to extract 
the two-quark magnetic moment operators results, if eq. I|3.30|l is cast in the form 



lim 



(3.31) 



which is similar to that obtained in ref. jlH]. By means of eq. (j3.3HI . it is now possible 
to consider the two-quark current operators for the scalar confining and vector OGE 
interactions, as calculated from the diagrams in Fig. 13.11 in paper V and ref. The 
two-quark current operator associated with the scalar confining interaction is then of the 
form 



Jox(g>fci,fc2) 



Qi-Pi 



(cTi + (T2) X q 
((Ti — 0-2) X q 





, Q*2 


2ml 


2m| 


Ql 


Q*2 


2ml 


2m| 



(3.32) 



where the variables Ql and Q2 are defined as Ql — Vc{k2)Qi and Q2 — Vc{ki)Q2, respec- 
tively. The corresponding current operator for the OGE interaction may be expressed 
as 



3jiiQ^ki,k2 



i<Ti X fca 2P2 + i(T2 X fc2 

2ml 



2mim2 



i(T2 X ki 

2m| 
2Pi + 



(3.33) 
icTi X ki 



2mim2 



with Ql — Vg{k2)Qi and Q2 = Vg{ki)Q2- As the above equation depends only on ki and 
^2, the OGE magnetic moment operator is most conveniently calculated using eq. (|3.31(l . 
The corresponding spin-flip operators for transitions between 5-wave quarkonium states 
have been obtained in paper V as 



Ql Q2 



((Ti - (T2) + 



Ql , Q2 



((Ti + (T2) 



(3.34) 



for the scalar confining interaction, and 



eVgir) 



Ql Q2 Ql — Q2 



mim2 

Ql , Q2 



{rri 

Qi- 



- (T2 

Q2 



m\m2 



{rri 



f^2, 



(3.35) 



for the OGE interaction. 
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For equal constituent quark masses, eqs. H3.34|l and (|3.35l) reduce to the expressions given 
in ref. 021 • Note that the presence of a spin-flip term in the OGE operator (|3.35() is solely 
a consequence of the difference in mass between the constituent quarks, and will thus not 
contribute to the Ml widths of the charmonium and bottomonium states. Similarly, the 
terms that are symmetric in the quark spins vanish for equal mass quarkonia. However, 
in the case of the system, these terms will contribute to the magnetic moment of the 
cb system. Also the spin-flip Ml transitions in the system will receive a contribution 
from the OGE operator. 



3.3 Widths for radiative decay 

The widths for El dominated transitions of the type Xcj J / ip j oi ip' ^ Xcj 7 have in 
paper VI been calculated according to 



r 



S 



q a 



9 



\Mo\ 



\M2 



(3.36) 



where J/ is the total angular momentum of the final quarkonium state, and q is the 
momentum of the emitted photon. The widths for ij}' XcJ 7 with J = 0, 1, 2 are then 
expected to scale as 1 : 3 : 5 respectively, but that result is highly modified by the large 
hyperfine splittings of the QQ with L = 1. The statistical factor Sfi is defined as in 
ref. |5()j and assumes the values Sfi — 1 for a triplet-triplet transition and Sfi = 3 for 
a singlet-singlet transition. On the other hand, in paper VI the widths for transitions 
between D- and P-wave states were obtained from 



^ „ 2J/-t-l 3 Mf .^^ 



S 



18 



J T) 



Jd 
1 



(3.37) 



where Jd and Jp are the total angular momenta of the D- and P-wave states, respec- 
tively. The values of Sfi are then given by the above Wigner 6j symbol. Note that the 
triangularity of the 6j symbol requires that \Jd — Jp\ = 1 or 0. Transitions that change 
the value of J by more than one unit are thus forbidden. In eqs. H3.36|l and (|3.37() . M.^ 
and A^2 denote radial matrix elements for S- and _D-wave photon emission, respectively. 
The radial matrix element for S-wave emission receives contributions not only from the 
impulse approximation, eq. H3.9|l . but also from the confinement and OGE operators in 
eq. (|3.16|l . That matrix element may thus be expressed as 



Mo = / drru/(r)Mj(r)jo (y) 



r dr 



(3.38) 



where Ui and Uf are the reduced radial wave functions for the initial and final heavy 
quarkonium states. Similarly, the matrix element for D-wave emission, which vanishes 
in the El approximation, is of the form 



M2 = iQ)mJ^ drrw/(r)u,(r) j2 (y) 



(3.39) 



The contribution from this matrix element is usually very small as the product gr <C 1 
for typical values of the photon momenta and quarkonium radii. Therefore, that matrix 
element has not been included in eq. H3.37|l . 
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The impulse approximation charge factor (Q)j^, and the exchange charge factors {Q)^ 
for the scalar confining interaction and {Q) g for the OGE interaction that appear in 
eqs. (I3.38|) and H3.39|l are of the form 



ia 



Qi 1 



m2 



%ra( J m,i + m2 



Q2 1 



8to| 



mi 



mi + m2 



(3.40) 



for the impulse approximation, where the quark charge operators have been multiplied 
with the Darwin- Foldy terms from eq. H3.10|l , and 



(Q>c = 



m2 



mi 



4mj^ mi + m2 4m2 mi + m2 
1 ^ 2 {Sf\rTi-a2\S,) 
mi 3 m2 



(3.41) 



Ami 



4m| 



1 

m2 



2 {Sf\cT^-a2\S-, 
3 



mi 



for the charge factors that are associated with the scalar confining and OGE interactions, 
respectively. In the spin dependent terms of eqs. (|3.41|l . Si and Sf denote the total spins of 
the initial and final quarkonium states. For triplet-triplet and singlet-singlet transitions, 
{Sf\(Ti ■ (T2\Si) = +1 and —3, respectively. The charge factor in eq. (|3.39|l is defined 

as ~ limg_^o {Q)i\- This is permissible since the Darwin- Foldy and exchange 

charge terms are very small compared to the dominant dipole contribution, which in 
itself is already insignificant because of the suppression by the j2 function in the matrix 
element. 



The expression for the width of a spin-flip Ml transition between S'-wave heavy quarko- 
nium states can be written in the form 



25, 



16 3 Mf .^^ . 



(3.42) 



where Ai.y denotes the radial matrix element for an Ml transition and Si is the total 
spin of the initial state. That matrix element consists of the relativistic impulse approx- 
imation, scalar confining and OGE components, according to 



M 



RIA 



M 



Conf 



(3.43) 



where the matrix element in the relativistic impulse approximation is defined according 
to 



M 



RIA 



drdr'rr'uf{r')ui{r) / dP P"^ 



1 



Qi 
mi " 



/7 



12 r-f 
^J2 



m2 



k {r'P)k (rP), 

(3.44) 

where the factors /7 are given by eq. (|3.28(l . The matrix elements associated with the 
scalar confining and vector OGE interactions, which have been shown to be large in 
papers I,V and VI are, in the nonrelativistic approximation, of the form 



A/fConf 



druf{r) Ui[r) 
druf{r) Ui(r) 



Vc{r) 



Qi Q2 



m. 



Q2 Qi — Q2 



mim2 



(3.45) 
(3.46) 



In particular, eq. H3.45|l is shown, in the next section, to provide an explanation for the 
experimental width of the Ml transition J — s- 77C7. 
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3.4 El and Ml transitions in heavy quarkonia 

A detailed comparison of the numerical results obtained in papers V and VI with exper- 
imental results and other theoretical calculations is instructive, as there are issues with 
several of the El and Ml transitions that are not readily apparent by casual inspection 
of the large amount of numerical data presented in those papers. This is even more 
important as the branching fractions for various transitions are typically better known 
than the total width of the initial state. With this in mind, the most important ones of 
the Ml and El transitions given in Tables mi and rOl are discussed below. 

3.4.1 The Ml transition J/ip rjc^ 

The major importance of this Ml transition, from both experimental and theoretical 
points of view, has been stated e.g. in the review of ref. 18'. The experimental width of 
1.14 ± 0.39 keV has been difficult to explain theoretically, since nonrelativistic calcula- 
tions overestimate this width by a factor ~ 3. A possible solution for this overprediction, 
which was already hinted at in ref. |47j , is presented in Table 13.11 where the exchange 
current contribution from the scalar confining interaction brings the width down to the 
desired level. The importance of such negative energy components for the transition 
J /ij) r]c"f has also been demonstrated within the framework of the instantaneous 
approximation to the Bethe-Salpeter equation in ref. '17' and within the Schrodinger 
approach in paper I. 

If the entire QQ potential had effective vector coupling structure, which has often been 
suggested in the literature j51| . then no exchange current contributions would arise, as 
a vector interaction contributes a significant spin-flip operator only if the quark and 
antiquark masses are unequal, and agreement with experiment would thus be excluded. 
Furthermore, a large family of effective vector confining interactions have been shown to 
be inconsistent with QCD by ref. [S^. However, it has also been shown in paper I that 
an expansion of the RIA spin-flip operator to order v'^/c^ overestimates the relativistic 
correction to the static (NRIA) result, which originally led to an opposite conclusion [Tfj 
concerning the usefulness of a scalar two-quark spin- flip operator. It was also suggested 
that the charm quark might possess a large anomalous magnetic moment, but that 
possibility has apparently not been substantiated. 

3.4.2 The Ml transition V' ^ ^c7 

This nonrelativistically forbidden Ml transition has also proved challenging to explain 
theoretically, since the (near) orthogonality of the quarkonium wave functions renders 
the results hypersensitive to small effects. In the recent calculation by ref. ^121, where 
good agreement with experiment was found for J/iJj — s- rjc 7, the width for ip' rjcj was 
however overpredicted by almost an order of magnitude. It is shown in Table mi that the 
Ml model employed in paper VI gives a width of ~ 1.1 keV for that transition, which is 
close to the upper uncertainty limit of the current empirical result 0.84 ± 0.24 keV |46| . 
That such a favorable result is obtained depends on several factors in the present work, 
such as the employment of ip ' and rjc wave functions that model the spin-spin interaction 
in the S'-wave. 
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The choice of approximation for the Ml matrix clement is also important in this respect. 
The amplitude (|3.17l) has the additional advantage of allowing the use of a realistic pho- 
ton momentum in the expression H3.42|l for the Ml width. Also, this treatment yields the 
same spin-flip operators as in the calculation of the exchange magnetic moment operators 
in ref. |49j . where the rigorous Ml approximation was used. Furthermore, the Ml approx- 
imation has been taken to affect the entire factor in brackets in eq. (|3.17|) . which leads 
to the elimination of the photon momentum q from the RIA matrix element H3.44|l . If 
the exponentials were separated from the current operators in eq. (|3.17() , then the width 
for ip' — > ?7c 7 would be overpredicted by a factor ~ 4. However, if spin-averaged wave 
functions were employed, as in paper I, then the conclusion would be exactly the oppo- 
site; In that case the present treatment would lead to unfavorable results. As seen from 
Table IXTl the exchange current operator associated with the scalar confining interaction 
gives the main contribution to the width for tp' ^ Tjc^ within this calculation. 

The present treatment of the Ml approximation may be regarded as consistent since 
it leads to the correct spin-flip operators and simultaneously allows a realistic photon 
momentum to be used. However, the large photon momentum introduces an additional 
uncertainty, which involves boosts on the QQ wavefunction in the flnal state, an effect 
which is yet to be considered. 

3.4.3 Other Ml transitions 

In principle, the width for T — ^ 7/h 7 could be predicted with much better accuracy than 
the corresponding one in the cc system, because of the large mass of the bottom quark. 
In particular, the exchange current contribution from the scalar confining interaction 
is much smaller than for cc. The largest uncertainty is introduced by the unknown 
photon momentum for the T — > ?7{, 7 transition, as the mass of the rjh state is not known 
empirically. As realistic models of the spin-spin splittings for S'-wave quarkonia give an 
rjb mass around 9400 MeV, then the width for T ^ y^h 7 is likely to be less than 10 eV, 
as given in Table IXTl 

In addition to the Ml transitions discussed above, predictions have also been given in 
paper VI for Ml transitions between QQ states below the thresholds for fragmentation. 
Among these is the transition ^ ' —^ r?^ 7, which is similar to J ftp — >■ Vc 7- As recent 
experimental results indicate that the mass of the 77^ is much higher than previously 
thought '33', then the amount of phase space available for ip' ^ rj'^j is also smaller. The 
predicted width for that transition is thus significantly smaller than the values suggested 
by previous work ^^l- The width for 77^ J/i^l is also sensitive to particulars of the 
model because of cancellations in the matrix element and the large photon momentum 
involved. The results in Table im suggest that the width for this transition should be 
around 2 keV. As the experimental situation concerning the rj'^ continues to improve, 
then the width for ry^ — > J/iJj 7 may possibly be measured in the near future. 

In the case of the bb system, the number of measurable Ml transitions is larger since the 
3S states of bottomonium lie below the threshold for BB fragmentation. These widths 
are difficult to predict and provide an important test for models of the Ml transitions. 
The results of paper VI suggest that the widths for transitions which do not change the 
principal quantum number of the quarkonium state should be highly suppressed, whereas 
the widths for transitions from excited rji, states to the T ground state are predicted to 
have larger widths of about 100 eV. 



Chapter 3. Electromagnetic Transitions 



30 



Table 3.1: The Ml transitions between low-lying S'-wave states in the charmonium (cc), 
bottomonium (bb) and {cb, cb) systems. Experimental data j46) is available only for 
J/%Ij ^ rjc^ and -0' — > '7c7- The quoted photon momenta have been obtained by 
combination of the empirical masses of the spin triplet states with splittings given by the 
Hamiltonian model of paper VI in Table ITTl Note that the OGE interaction contributes 
only to Ml transitions between Be states. 



Transition 






Matrix element [fm] 








Width 




NRIA 




ria 


Exch 




NRIA 


RIA 


RIA+Exch 


J/ip ^ ric-y 
: 116 McV 


4.356 


10" 


-2 


3.762 


10-2 


-8.724 


10" 


3 


2.85 


2.12 


1.25 keV 

cx: 1.14 ± 0.39 


: 639 McV 


3.985 


10- 


-3 


-5.14 


10-* 


2.826 


10" 


3 


3.35 


0.06 


1.13 keV 

cx: 0.84 it 0.24 


V" ' V'cl 
: 46 MoV 


4.344 


10" 


-2 


3.735 


10-2 


-1.870 


10" 


2 


0.18 


0.13 


0.03 keV 


q^ : 502 McV 


-4.271 


10" 


-3 


-7.584 


10-3 


5.206 


10" 


3 


5.89 


18.6 


1.83 keV 


T %7 
q^ : 59 McV 


-6.71 


10" 


-3 


-6.39 


10-3 


2.41 


10" 


4 


9.2 


8.3 


7.7 eV 


q-, : 603 McV 


-3.94 


10" 


-4 


-1.44 


10-* 


-8.76 


10" 


5 


31.8 


4.3 


11.0 eV 


T' ^ < 7 
q^ : 25 McV 


-6.70 


10" 


'3 


-6.39 


10-3 


5.45 


10" 


4 


0.70 


0.64 


0.53 eV 


'?{, ^ Tf" 7 
q-, : 530 McV 


4.18 


10" 


-4 


6.30 


10-* 


-1.31 


10" 


4 


71.5 


162 


102 eV 


T" ^ r,'^ 7 
q-, : 16 McV 


-6.70 


10" 


'3 


-6.35 


10-3 


8.02 


10" 


4 


0.18 


0.16 


0.13 eV 


T" ^ri'^j 
q-, : 350 McV 


-3.59 


10" 


-4 


-1.11 


10-* 


-1.55 


10" 


4 


5.3 


0.5 


2.9 eV 


r" ^Vb-y 

q^ : 910 McV 


-2.10 


10" 


-4 


-6.59 


10-5 


-3.77 


10" 


5 


30.2 


3.0 


7.3 eV 


< ^T'7 
q^ : 311 McV 


3.96 


10" 


-4 


6.05 


10-* 


-2.25 


10" 


4 


13.7 


32.0 


12.6 eV 


< - T7 
q^ : 842 McV 


2.05 


10" 


-4 


3.02 


10-* 


-4.68 


10" 


5 


68.7 


149 


106 eV 


B* ^ Bc7 

q-, : 53 McV 


1.851 


10" 


-2 


1.496 


10-2 


0.326 


10" 


3 


50.0 


32.6 


34.0 eV 


B*' -> Bc7 

q^ : 576 McV 


1.015 


10" 


-3 


-1.437 


10-3 


2.524 


10" 


3 


179 


360 


206 eV 


B*' - 5^7 
q-, : 22 McV 


1.849 


10" 


-2 


1.480 


10-2 


-5.186 


10" 


3 


3.61 


2.31 


0.98 eV 


B', ^ B* 7 

q^ : 507 McV 


-1.067 


10" 


-3 


-3.089 


10-3 


1.028 


10" 


2 


411 


3440 


39.5 eV 
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3.4.4 The El transitions Xcj ^ Jj'^^ and ^' XcJl 

The El transitions from the spin-triplet P-wave states are in principle the simplest to 
predict accurately, as the wave functions involved do not contain any nodes. Although 
the empirical data from ref. 32 has suggested that the El widths are generally over- 
predicted |44| , that discrepancy is apparently resolved by the new data presented in the 
2002 edition of the PDG 46 . However, the results in Tablc rT^ indicate that the rigorous 
El approximation, with q — Mi — Alf, overpredicts the widths by a factor ^ 2. If the 
El approximation is removed, the recoil of the J/^ can be accounted for, in which case 
that overprediction is eliminated. 

Prediction of the widths for 7 transitions from the ijj' state has proved to be difficult, 
as the El approximation typically overpredicts the widths by at least a factor ^ 2. The 
present empirical data on the ip' suggests that the widths for i/i' — > XcJ7 should 
be around 25 keV. As demonstrated in paper VI, the El approximation yields widths 
in excess of 40 keV. This is puzzling, since recoil effects are small and cannot explain 
this overprediction. Also, it is seen by inspection of Table that the predicted relative 
widths also do not agree well with experiment, although the experimental uncertainties 
are considerable. Not surprisingly, the matrix elements in Table 13.21 reveal that these 
transitions are very sensitive to small hyperfine effects in the QQ wave functions, which 
has also been demonstrated in ref. It is therefore conceivable that small modifica- 
tions of the QQ wave functions may be sufficient to solve this overprediction. It should 
also be noted that significant reductions of the El widths were achieved in ref. by 
consideration of closed cq — qc fragmentation channels. 

3.4.5 The El transitions from the Xbj states 

The calculated widths for the Xbj — >■ ^7 transitions agree rather well with those of the 
other models presented in Table ITSI although they appear to be slightly larger. If the 
calculated El widths are used to predict the total widths of the Xbj states, then it is 
found that the width of the Xb2 should be 164 ± 22 keV and that of the Xfci about 93 ± 22 
keV. Similarly, the calculated El width of the Xfeo suggests that the total width of that 
state is at least ^ 440 keV. This situation is similar to that observed for cc 32_, where 
the Xc2 is wider than the Xci by about a factor 2. 

The El transitions from the Xf,j(2P) states in bottomonium provide a useful test for 
theoretical models since experimental data is available on all six branching fractions |32| . 
even though the total widths of the Xfcj(2P) states are not known. These data indicate 
that the widths for transitions to the T should be about one half of those for transitions 
to the T(2S'), even though much more phase space is available for the former. Indeed, 
it can be seen from Table rOl that spin-averaged wave functions do not provide a good 
description of the experimental branching fractions even though the hyperfine splittings of 
the Xbj{'2P) states are small. On the other hand, much better agreement with experiment 
is obtained if the hyperfine effects are accounted for by the QQ wave functions. The 
calculated widths for X6j(2P) T{2S)j may be used to estimate the total widths 
of the X6,7(2P) states from the known branching fractions. The predicted width of the 
Xb2 (2P) state is then 100 ± 15 keV, while that of the Xbi (2P) is 72 ± 14 keV. The Xbo (2P) 
state appears to be significantly broader, but because of the large errors in the reported 
El branching fractions, only a rough estimate of 267 ± 140 keV is possible. 
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Table 3.2: The El dominated transitions between low-lying states in the charmonium (cc) 
system, together with the empirical data given by ref. |B2j. The column "lA" contains the 
matrix element H3.38|l in the impulse approximation, while in the column labeled " DYN" , 
the exchange charge contributions have been included. Note that a good experimental 
candidate |S2 for the ^Di state is the V (3770) resonance. 



Transition 




Mo [fm] 




M2 


fm] 




Width 




lA 


DYN 


El 








El 


DYN 


Xc2 J/V'7 
: 429 MeV 


0.2389 


0.2442 


0.2632 


7.145 


■ 10~ 


3 


558 keV 


343 keV 

cx; 389 ± 60 


Xcl * 'J 1 ^^1 

■ 390 MeV 


0.2464 


0.2519 


0.2673 


5.729 


■ 10" 


3 


422 keV 


276 keV 

cx: 290 ± 60 


XcO J 1^1 
q^ ■ 303 McV 


0.2556 


0.2612 


0.2701 


3.345 


■ 10" 


3 


196 keV 


144 keV 

exp: 165 ± 40 


i)' XcOl 
q^ ■ 261 McV 


-0.2685 


-0.2686 


-0.2840 


-6.106 


■ 10- 


-3 


44.6 keV 


33.1 keV 

cx: 26.1 ± 4.5 


V' ' ^ Xcl 7 
q^ ■ 171 McV 


-0.3126 


-0.3126 


-0.3202 


-3.028 


■ 10- 


-3 


45.8 keV 


38.7 keV 

cx: 25.2 ± 4.5 


'Z' ' ^ Xc2 7 
q^ : 127 McV 


-0.3440 


-0.3442 


-0.3489 


-1.871 


■ 10- 


-3 


37.1 keV 


33.1 keV 

ox: 20.4 ± 4.0 


he ^ rjc'J 
q^i ■ 493 McV 


0.2098 


0.2091 


0.2289 


7.377 


■ 10- 


3 


661 keV 


370 keV 


^ he 7 
q-y : 125 McV 


-0.3420 


-0.3424 


-0.3465 


-1.618 


■ 10- 


-3 


61.5 kcV 


55.0 keV 


i)" XcO 7 
qy : 577 McV 


-0.0456 


-0.0450 


-0.0199 


0.926 


■ 10" 


2 


2.69 keV 


9.86 keV 


V'" Xcl 7 
q^ : 494 McV 


-0.0306 


-0.0298 


-0.0033 


1.016 


■ 10- 


2 


0.13 kcV 


9.57 keV 


^" Xc2 7 
q^ : 455 MeV 


-0.0168 


-0.0161 


0.0123 


1.099 


■ 10- 


2 


2.38 keV 


5.75 keV 


r^x'coJ 

q^ : 153 McV 


-0.4315 


-0.4315 


-0.4497 


-7.344 


■ 10- 


-3 


21.3 kcV 


17.8 keV 


r ^ X',1 7 
q^ : 125 McV 


-0.4860 


-0.4861 


-0.4995 


-5.399 


■ 10- 


-3 


42.6 keV 


37.3 keV 


r^x',2^ 

q^ : 109 MeV 


-0.5280 


-0.5283 


-0.5391 


-4.367 


■ 10- 


-3 


53.7 keV 


48.2 keV 


^Ds Xc2 7 
q^ : 227 McV 


0.4164 


0.4194 


0.4353 








243 keV 


192 keV 


^D2 Xc2 7 
q^ : 221 MeV 


0.4188 


0.4219 


0.4367 








56.5 keV 


45.2 keV 


Xcl 7 
q^ : 263 McV 


0.3920 


0.3953 


0.4145 








262 keV 


198 keV 


^Di Xc2 7 
q^, : 206 MeV 


0.4216 


0.4246 


0.4372 








5.06 kcV 


4.13 keV 


^Di Xcl 7 
q^ : 248 McV 


0.3963 


0.3997 


0.4164 








123 keV 


94.9 keV 


^Ol XcO 7 
9.^ : 336 McV 


0.3578 


0.3619 


0.3889 








370 keV 


251 keV 
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3.4.6 The El transitions from the T states 

The experimental situation concerning the T(2S) XbJl transitions has lately become 
more uncertain since the total width of the T{2S) as reported by ref. (32,, originally 
given as ~ 27 keV, has increased over time and now stands at 44 ± 7 keV. This situation 
is analogous to that for the ip\ which has undergone a similar, albeit smaller, increase. 
This has made the model predictions in Table rOl which originally fitted the experimental 
data very well, much less satisfactory. It is therefore very difficult to judge the quality of 
any given prediction until the experimental situation is stabilized. Still, it is noteworthy 
that the calculation of paper VI gives slightly better agreement with experiment than 
the other models in Table 

As the reported total width of the T(3S') state [22, 26.3 ± 3.5 keV, is better known 
than that of the T(25') state, then it is expected that systematic uncertainties in the 
reported experimental results for T(3S') Xbj{'^P)l should be smaller than for the 
analogous T(2S') XbJl transitions. By inspection of Table it can be seen that 
the T(3S') Xbj{'^P)l transitions are generally rather well described by a number of 
models, although the calculation of ref. j47j . where spin-averaged wave functions were 
employed, apparently underpredicts the empirical widths. Also, the results of ref. |54| 
appear to compare slightly more favorably with experiment than those of the calculation 
in paper VI. 

While the T(3S') Xb.ii'^P) 7 transitions are relatively well described by different mod- 
els, the situation concerning the T(3S') Xbj 7 transitions remains unsettled because of 
a strong cancellation in the El matrix element. However, experimental detection of these 
transitions may be a formidable task since the widths are an order of magnitude smaller 
than those of any previously measured El transition in the bb system. Inspection of the 
results in paper VI reveals that within the dynamical model, the width for T(3S') — > Xbo 7 
should be the largest and that for T(3S') — > Xb2l the smallest. It is encouraging that 
the same pattern is also predicted in Table 13.21 for the analogous transitions in the cc 
system, where the widths are much larger relative to the other El transitions. 

3.4.7 Other El transitions 

The widths for El transitions from the ip{3S) state have also been given in Table POl 
Those results suggest that the transitions to the (2P) states should have widths that are 
comparable to those for the ip' — > XcJ7 transitions. On the other hand, the widths for 
the ipi-^S) — ^ XcJ 7 transitions are predicted to be smaller by factors 3 — 4. The empirical 
detection of any of these transitions will probably be difficult since the ?A(3S') state lies 
above the DD threshold. However, since the spin singlet he state is well below threshold, 
then the photon produced in the he — > rjc"/ transition may be detected in the near future. 
The dynamical model yields a width of 370 keV for this transition, which is then the 
largest El width in the cc system, although the difference between the El approximation 
and the dynamical model is large for that transition. 

The El transitions from the "^Di state are also of particular interest, as that state proba- 
bly corresponds to the empirical -0(3770) resonance. The predictions of paper VI suggest 
that the transitions to the Xci and Xco states should be detectable by experiment, whereas 
that to the Xc2 state is highly suppressed by the statistical factor Sfi. 
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Table 3.3: Comparison of the predicted El widths in the bottomonium (bb) system with 
those of other models that use a scalar confining interaction. All widths are given in 
keV. The experimental widths have been extracted from the branching fractions and 
total widths reported by ref. 132! . 





GS (ref. 1171) 


GZ (ref. [511) 


paper VI 


Exp (ref. ED) 


Xb2 - 


T7 


33.0 


33.8 


36.0 


22 ± 3% 


Xbl ^ 




29.8 


30.4 


32.5 


35 ± 8% 


XbQ - 


T7 


25.7 


25.3 


26.6 


< 6 % 


r^/ 


> XbO 7 


0.73 


0.76 


1.01 


1.7 ±0.5 keV 




* Xbl 7 


1.62 


1.37 


1.80 


3.0 ± 0.7 keV 


X' — 


* X&2 7 


1.84 


1.45 


2.03 


3.1 ± 0.7 keV 


ylo - 


^Vj 


12.9 


16.2 


16.4 


16.4 ±2.4 % 


Xbl - 


T'7 


11.9 


14.7 


15.1 


21 ± 4 % 


X'bo ' 


T'7 


10.6 


12.3 


12.3 


4.6 ± 2.1 % 


Xb2 - 


-.T7 


18.2 


10.4 


11.4 


7.1 ± 1.0 % 


Xbl - 


-.T7 


11.8 


7.51 


8.40 


8.5 ± 1.3 % 


Xm - 


-.T7 


6.50 


3.57 


4.93 


0.9 ± 0.6 % 


Y" — 


^ Xho7 


0.114 


0.029 


0.15 




Y" — 


^ X617 


0.003 


0.095 


0.11 




Y" — 


^ X62 7 


0.194 


0.248 


0.04 




_ 


^ Xho7 


1.09 


1.30 


1.14 


1.4 ±0.3 keV 




Xfci7 


2.15 


2.34 


2.12 


3.0 ±0.5 keV 


Y" „ 


^ Xfc2 7 


2.29 


2.71 


2.50 


3.0 ± 0.6 keV 



The determination of the photon momenta for transitions in the bottom-charm sys- 
tem has to rely on model predictions for the masses of the cb states. The uncertainty 
introduced by this is, however, rather small for the El transitions, as the model predic- 
tions for the major level splittings agree with each other to a large extent Inspec- 
tion of the results in paper VI reveals that the predicted widths are similar to those 
obtained by ref. [SUI, although differences exist for transitions like B*{2S) S*q7 and 
5*2 (2P) B*{2S) 7, where the widths are sensitive to the effects of the hyperfine com- 
ponents of the QQ interaction. It is noteworthy that while the predicted widths for the 
B*j -^5*7 transitions agree rather well with those from ref. [HOI , there is a signifi- 
cant disagreement for Bd — > Be 7. When the somewhat different photon momenta are 
accounted for, this disagreement amounts to about a factor ~ 3. 

An issue not considered in paper VI is the spin mixing of the L — 1 states with J = 1, 
that is due to the antisymmetric spin-orbit interaction which was not included in the QQ 
interaction Hamiltonian. This mixing, which was considered in ref. |50| . has the effect 
of allowing "spin- flip" El transitions of the type B*^ — » i?c7- However, the widths for 
such "forbidden" transitions were found in ref. [101 to be typically suppressed by a factor 
~ 100 relative to the "allowed" ones considered in paper VI. 
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3.5 Ml transitions in heavy-light mesons 



The Ml widths of the spin-flip Ml transition between the vector and pseudoscalar states 
in the charm mesons have been calculated in paper V, where the Qq interaction was 
modeled as a scalar confining + OGE potential. A comparison of the results presented in 
Table im is instructive since the total width of the D^* state has recently been measured 
by the CLEO collaboration [SSj. The width is reported as T{D^*) = 96±4±22 keV, where 
the latter error represents the systematic uncertainty. The reported |32j branching ratio 
of 1.6 ±0.4% for radiative decay then gives a width of 1.5 ±0.6 keV for the Ml transition 
D^* D^7. Here most of the uncertainty can be traced to the systematic errors of the 
experimental result. Such a comparison shows that the results of Table reproduce 
the empirical width fairly well for a range of values of the light constituent quark mass 
niq. The value = 450 MeV corresponds to the potential model of ref. 38 , while the 
value ruq = 420 MeV has been suggested in ref. |^ . Indeed, for a light constituent quark 
mass of 420 MeV, a width for radiative Ml decay which is close to 1.5 keV is reproduced. 
Values close to that are also favored by the analysis within the framework of the Gross 
equation by ref. |56| . 





NRIA 


RIA 


RIA + Exch 


450 MeV 


21.1 keV 


8.86 keV 


8.95 keV 


420 MeV 


23.5 keV 


9.18 keV 


9.89 keV 


390 MeV 


26.4 keV 


9.52 keV 


11.1 keV 





NRIA 


RIA 


RIA + Exch 


450 MeV 


0.58 keV 


9.4 


10-3 


keV 


1.09 keV 


420 MeV 


0.79 keV 


1.5 


10-2 


keV 


1.43 keV 


390 MeV 


1.07 keV 


2.2 


10-2 


keV 


1.90 keV 




NRIA 




RIA 




RIA + Exch 


560 MeV 


0.18 keV 


2.6 


10-4 


keV 


0.38 keV 


530 MeV 


0.26 keV 


3.9 


10-5 


keV 


0.49 keV 


500 MeV 


0.36 keV 


8.6 


10-4 


keV 


0.64 keV 



Table 3.4: Numerical results 
from paper V for the Ml transi- 
tions between ground state vec- 
tor and pseudoscalar mesons in 
the D and Dg systems, for 
different values of the light 
and strange constituent quark 
masses. The charm quark mass 
of 1580 MeV corresponds to the 
potential model of ref. In 
the right-hand column, the two- 
quark exchange current contri- 
butions from the scalar confin- 
ing and OGE interactions have 
been added to the RIA result. 



Even though the total width of the neutral Z)"* meson has not yet been determined |32| . 
considerable information about the expected width for D"* — > D°7 may be extracted 
from the reported branching fraction of 38.1 ± 2.9% since the corresponding width 
for pion emission can be constrained by means of the empirically determined width of the 
D^* and model calculations of the pionic transitions in D mesons 'SH', 'Ff. If one notes 
that the branching fraction of D^* D^n^ is reported as 67.7±0.5% which implies 
a width for this transition of ~ 65 ± 14 keV, then it is found from the model calculation 
of paper III that this corresponds to a width of ~ 40 ± 10 keV for D°* D°Tr°. As the 
relative branching fractions for 7r° and 7 emission by the D'^* are well known |221, the 



Chapter 3. Electromagnetic Transitions 



36 



best estimate for the width of I?"* D^^ is ^ 25 keV, which is close to that preferred 
by ref. |56| . There remains, however, a considerable uncertainty of ^ ±10 keV from the 
systematic errors in the empirical measurement of T[D^*). 

The results of paper V in Table underpredict this expectation by about a factor ~ 2, 
which underlines a basic weakness of the present approach to the Ml transitions in the 
D mesons. Firstly, as the two-quark spin-flip operators of eqs. (|3.34() and H3.35|l are given 
for the nonrelativistic approximation, a relativistic treatment of those operators will, in 
general, lead to a weakening of the two-quark contributions to the matrix element for an 
Ml transition. It is thus entirely possible that a fully relativistic treatment will render 
the two-quark contributions too weak to account for the experimental data on Ml decay 
of the meson as well. This conclusion is in line with that reached in ref. [SB|, where 
the discrepancy between the model and experiment was parameterized in terms of a large 
anomalous magnetic moment for the light quark. 

Another interesting possibility discussed in paper V suggests itself, in view of the prob- 
lems of fitting the spectra of the Qq mesons within the framework of the BSLT equation 
using a OGE interaction alone, namely that the instanton induced interaction for Qq 
systems given in ref. |3(J| may also contribute a significant two-quark current. The inter- 
action of ref. |3U| . which was found to be short-ranged, attractive and with negative sign, 
has scalar coupling to the light constituent quark. It has been noted in paper V that 
such an interaction adds up constructively with the OGE contribution, and counteracts 
that from the scalar confining interaction. An overall favorable effect on the widths for 
Ml transitions may thus be obtained, which may be inferred from the matrix elements 
given in paper V. 



Chapter 4 

Pionic Transitions 



The pionic transitions in the heavy-hght D mesons are instructive, as they can provide 
direct information on the strength of the couphng between pions and hght constituent 
quarks. Furthermore, as the charm quark in the D mesons does not couple to pions, the 
decay mechanism is governed by the pion couphng to the hght constitutent quark alone. 
As a first approximation, the pion-light constitutent quark coupling can be taken as 
being independent of the quark-antiquark interaction in the D meson. It has been shown 
in paper II, where the pion emission was modeled in terms of the chiral pseudovector 
Lagrangian, that while this assumption is reasonable for the axial current term, it leads 
to a large overestimate of the axial charge contribution. In this case two-quark pair 
terms analogous to those that are required for a realistic description of the Ml transition 
J/tp ^ rjcl may reduce the axial charge amplitude to a realistic level. An analogous 
suppression of the S'-wave pion transitions was achieved in ref. ,56i within the framework 
of the Gross quasipotential reduction. 

The empirical information on the widths and branching fractions of excited charm mesons 
is still very limited. Absolute values, albeit with large uncertainties, are known for the 
£'i(2420) and £'2(2460) mesons, and recently the width of the D*(2010)^, for which 
only an upper bound of 0.131 MeV was available earlier, has now been reported as 
96 ± 4 ± 22 keV by the CLEO collaboration [SSj . This result is shown to be consistent 
with values of the pion-quark axial coupling (7^ that are slightly smaller than 1. The 
observed widths of the L = 1 Z3i(2420) and D2(2460) mesons are also shown to be fairly 
well described, although a slight underprediction is expected as it has been demonstrated 
in paper III that two-pion tttt decay may also contribute significantly to the total widths 
of those mesons. In particular, the analogy with the i^|(1430) strange meson suggests 
that TTTT transitions may account for a significant fraction of the observed total widths. 

The flavor symmetry breaking decay mode D* — > Z?£,7r°, which has been considered in 
paper V is mainly due to a small isoscalar 77 meson component in the physical 7r° meson, 
as only the rj meson can couple to the strange quark in the Ds meson. The known ratio 
of the branching fractions for D* Dgir'^ and D* Dgj may be used to extract the 
coupling of 77 mesons to strange quarks, once the value of the n'^—r] mixing angle is known. 
Applied to the quark model for the baryons, an rjNN pseudovector coupling constant, 
small enough to be consistent with the phenomenological analysis of photoproduction of 
the 77 on the nucleon and the reaction pp — > ppr], has been obtained in paper V. 
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4.1 The amplitude for pion emission 



In paper II, the emission of pions from a D meson was described by the pseudovector La- 
grangian, which constitutes the lowest order chiral couphng for pions to hght constituent 
quarks: 



qq-n 



'(Pql5l^idf,T:aTatpq. 



(4.1) 



Here g\ denotes the axial coupling constant of pions to light constituent quarks, and 
is the pion decay constant, the empirical value of which is 93 MeV. The axial coupling 
constant is conventionally taken to be equal to, or somewhat less than, unity 58 . The 
Lagrangian H4.1|l yields the following transition amplitude for pion emission: 



g1 



(4.2) 



where ^ is an isospin factor, the value of which is for tt^ and 1 for 7r° emission. 
Decomposition of the above amplitude into axial current and axial charge components 
yields the pion emission amplitude in the impulse approximation. 



^'^2U V ^E'E„ 



9\ 



2m„ + E„ 



p2 _g2/4 



i{E' +mq){Eq+mq) 



jAEqE'^{Eq+mq){E'^ + mq) 



<Tq ■ q 



(4.3) 



where the first (axial current) term gives rise to the P-wave transitions D* Dtt and 
the D-wave transitions from the L = I charm mesons. On the other hand, the axial 
charge term leads to S'-wavc transitions from states with L = I. 




Figure 4.1: Irreducible two-quark contributions associated with the Qq interaction to the 
axial current and axial charge operators, with four-momentum variables defined as for 
Fig. 13.11 The pion always interacts with the light constituent quark, as the charm quark 
does not couple to pions. 
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As demonstrated in paper II, two-quark axial exchange current operators illustrated by 
Fig. 14.11 give large contributions to several tt transitions, in particular to those which 
involve the axial charge amplitude. In that case the axial exchange charge operator is of 
equal magnitude as the single quark operator, while the axial exchange current operators 
typically represent ~ 10% corrections to the impulse approximation result. If, in the 
static approximation, the axial current — (Aa,iAoa) of the light constituent quark 
is expressed as 

Aa = -gl(TqTa, (4.4) 

then the contribution to the axial current from a scalar confining interaction is, to lowest 
order in v/c, of the form (59) 



(4.5) 



where Vc{k2) is the Fourier transform of the scalar confining interaction. The above 
expression does not include the corrections from the canonical boost factors on the single 
quark spinors that are included in the single quark operator, eq. 1)4.3(1 . Moreover, a factor 
in the axial exchange current operator 1)4. 5|l arises as the static approximation to 
the propagator of the intermediate negative energy quark. 

Hence a more realistic evaluation requires that those factors are taken into account. 
For simplicity, the same spinor factors as for the single quark operator were used in 
paper II. The so obtained results indicate that the static approximation implies a very 
large overestimate of the axial exchange current contribution. Nevertheless, as shown 
in Tables 14.11 and 14.21 it serves to increase the calculated widths for pion emission. A 
complete calculation would also require consideration of the axial exchange currents 
associated with the short-range OGE or instanton induced interactions, which was not 
attempted in paper II. 

Both the scalar confining and OGE interactions contribute a two-quark operator to the 
axial charge amplitude in eq. 1)4.3(1 . These operators have been calculated in ref . [^01 and 
are of the form 

A'ir' = ^V,{k^)<Tq-PTa, (4.6) 



for the scalar confining interaction, and 



'^9 • + 9*^9 X (Tq • fc2 



(4.7) 



for the OGE interaction, which, when compared with the single quark axial charge 

^Oa = -— Tg-Pr, (4.8) 

reveals that the scalar confining interaction will tend to cancel out the axial charge 
component of the amplitude for pion emission. However, the calculations in paper II 
show that the OGE contribution is also large, although formally suppressed by a factor 
mq/MQ. It has also been shown in paper II that a relativistic treatment of the axial 
exchange charge operators will weaken them significantly. Although qualitative, these 
conclusions are in line with the experimentally small width of the L = I, J = I Di 
resonance. 
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4.2 The pionic widths of the D mesons 

The pionic decays of the D* mesons presented in Table 14.11 are intriguing since the 
emitted pions are extremely soft. Due to the very small phase space, the transition 
jj*o _^ D^TT^ is kinematically forbidden. The results appear to favor the value of 
= 0.87, although that is accidental since only an upper limit on the width of the D^* 
was known |32| when paper II was published. For this value of the total widths of 
the L — 1 D mesons in Table ^21 appear to be underpredicted, although tttt transitions 
may contribute significantly to these, as proposed in paper III. However, the empirical 
widths are still rather uncertain |46| , and have decreased over time. 

The results for the pionic widths of the excited D mesons are rather similar to those of 
ref. especially for the transitions D* — > Dtt, even though the Gross framework was 
employed in that paper. That calculation was restricted to the transitions allowed by the 
lowest order selection rules suggested by Heavy Quark Symmetry (HQS) jS^j whereas 
the present work uses the LS'-coupling scheme, which is more appropriate for equal mass 
quarkonia. The connection between the present calculation and the heavy quark limit 
remains as yet unexplored. 

Table 4.1: The calculated and experimental |46[ I55| pionic widths in McV for the D* 
mesons, corresponding to g\ = 0.87. The single quark approximation, with relativistic 
corrections is denoted RIA, and the result obtained upon addition of the axial exchange 
current contribution is denoted RIA + EXCH. The net results are also shown for = 1. 



Transition 


TT mom. 


RIA 


RIA + EXCH 


5^ = 1 


Experiment 




38.3 keV 


0.026 


0.029 


0.038 


0.029 ± 0.008 MeV 




39.6 keV 


0.056 


0.064 


0.084 


0.065 ±0.017 MeV 




43.1 keV 


0.036 


0.041 


0.054 


< 1.3 MeV 



Table 4.2: Calculated and empirical pion decay widths of the Di and mesons driven 
by the axial current and charge operators respectively, for g\ = 0.87. The empirical 
values are total widths 46 , which should mainly be due to pionic transitions to the 
ground state. The numbers in parentheses are the widths obtained without the axial 
exchange current contribution. The calculated values are also shown for g\ — l. 



Transition 


Current (RIA) 


Charge 


Total 


5^ = 1 


Experiment 


Di D*n 
Dtt 

Dl -> D*tj: 
Dt: + D*tt 


4.2 (3.4) 
8.1 (6.7) 
3.9 (3.1) 
11.9 (9.9) 


6.1 


10.3 MeV 
8.1 MeV 
3.9 MeV 
11.9 MeV 


13.6 
10.6 
5.1 
15.7 


I8.9tt l MeV 
25tr MeV 
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4.3 7T^ and 7 transitions from the meson 

Consideration of the coupling of the octet of light pseudoscalar mesons to the light {u,d,s) 
quarks yields, analogously to eq. (|4.1(l the couplings 



g1 



(4.9) 



For the pions and the t] meson, the empirical decay constants are /^r = 93 MeV and 
fri = 112 MeV, respectively, so at least for the decay constants SU{3) flavor symmetry 
is broken only at the 10% level. Combination of the chiral coupling (|4.9|l with the 
representation 



iPaXa = V2 





^+ 


K+ 






K" 








V 3 ' 






(4.10) 



gives the following definitions for the quark- level pseudovector coupling constants /i^gg, 



r ™7r q 



frjqq — 



9 A' 



fr]ss 



9 a- 



(4.11) 



The above relations then suggest that the magnitude of the coupling of rj mesons to 
M, d quarks should be one-half that of the rj coupling to strange quarks, independently 
of the 77 meson mass. In the static quark model the meson-quark coupling constants of 
eq. H4.11|l are related to the meson-nucleon coupling constants as 



NN 



■ fnqq, frjNN — fr] 



(4.12) 



Application of the relations (|4.11() together with eq. H4.9|l yields a coupling of 77 mesons 
to strange quarks in terms of /^s^ . The coupling of the tt*^ meson to the strange quark 
may thus be expressed in terms of the "effective" mixing angle 9m, which corresponds to 
the sum of the iT^—rj and contributions. The width for the process D* DsTt'^ can 

then be conveniently obtained from the expression for D^* — > D^tt^ given in paper II 
by the replacement g\/2f^ ^ fr^ssdm/m^, giving 



1 Md^ 



(4.13) 



if the tt'^ emission takes place at the strange quark. Here is a radial matrix element for 
pion emission. On the other hand, the width for the radiative Ml transition D* Dg^ 
was obtained, in paper V, as 



16 Md, 3 I , . ,2 



(4.14) 



where is a radial matrix element for Ml decay. By means of eqs. (|4.13f) and (|4.14|) . 
the ratio of the tt*' and 7 widths of the D* meson is then obtained as 



8 f2 /92 

o J r/ss rn 

97T m^^a 



\M,\ 



(4.15) 



where the dimension of |A^-y| is [MeV] 
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4.4 Estimation of fjjNN 

Through use of the empirical ratio of pion and photon momenta known to be ap- 
proximately 139/48 and the i] meson mass of 547 MeV one may solve for the coupling 
constant f,^ss to get 

As the ratio of the tt*^ and 7 decay rates is experimentally known, albeit with quite large 
errors, to be 0.062 ± 0.028 it is, given the rather well known value of 6m, possible 
to obtain an estimate for the coupling constant friss- However, if the charm quark also 
couples to 7r°, which is suggested by the empirically detected transitions ip' J/fpV 
and ip' — > J/t/i tt" , then eq. 14.16|) should be modified. The appropriate modifications 
are given in paper V, where it was found that a significant coupling of the tt*' to the 
charm quark will increase the value of fnNN- For that calculation, a matrix element for 
7r° emission by the charm quark is also required. 

The matrix elements required for eq. (|4.16|) were in paper V found to be 

= -1.22 • 10-2 fm 

for the Ml transition D* — > Ds 7, and 

A^^ = 0.794, Ml = 0.949, 

for 7r° emission by the strange and charm quarks. Assuming that in the 7r° emission by 
the 13* , the pion couples mostly to the strange constituent quark, eq. (|4.16() may be used 
directly together with the above matrix elements for the vr'^ and 7 transitions. Insertion 
of those matrix elements for a value of the 7r° — 77 mixing angle of 9m ^ 0.012 yields 
|/i)ss| ^ 0.70. If the uncertainties in the mixing angle and the empirical widths for the 
7r° and 7 transitions are taken into account, then the best estimate of paper V for the 
coupling of the 77 meson to strange constituent quarks is 

Ji-jss 'J- 1 — 0.3' 

In the above result, the negative sign is suggested by the relations in eq. H4.11|l . The static 
quark model then implies, through eq. H4.12|l . that the magnitude of the corresponding 
pseudovector 77-nucleon coupling constant fr^NN should be one half of this value. Thus 
one obtains the following final result for the 77-nucleon coupling: 



f-qNN — 0.35 _q;25 



This result should be compared with the value for fnNN or the equivalent pseudoscalar 
coupling constant q^nn = i'2'mN/'mri)friNN, which has been determined by phenomeno- 
logical model fits to photoproduction of 77 mesons on the nucleon [H^] . The latter value 
for fnNN is ~ 0.64. That value has also been found to be realistic in calculations of the 
cross section for pp pprj near threshold • Although the result obtained above for 
fr/NN has quite large uncertainties which are mostly of empirical origin, it still appears 
to be significantly smaller. A larger value for fnNN could, of course, be obtained by 
decreasing the mixing angle Om- 
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Transitions 



As the orbitally excited L ~ 1 £'i(2420) and 152(2460) charm meson states he weh above 
the threshold, not only for single pion but also for two pion decay, then it is likely that a 
significant fraction of their total widths are made up by tttt transitions to the ground state 
D and D* mesons. It is thus particularly instructive to obtain theoretical predictions and 
empirical information on the branching ratios for the latter decay modes. At present, 
however, the total widths of the Di(2420) and i?2(2460) states are known only within 
a very wide uncertainty range, and the remaining two members of the L — 1 multiplet 
have not yet been discovered. Paper III reports a calculation of the tttt decay widths of 
the excited L = 1 charm meson states, by extending a similar calculation of the widths 
of their single pion transitions in paper II. 

Two pion emission from radially excited heavy quarkonium (QQ) states empirically con- 
stitutes a significant fraction of their total decay widths Indeed, in the case of the 
ip' (or ■0(25')), the branching ratio is empirically as large as ~ 50%. As the charm quarks 
themselves do not couple to pions, the tttt coupling to heavy flavor mesons (or quarks) 
involves at least two gluons if not a glueball. A number of different theoretical approaches 
for the coupling of two-pions to heavy mesons have been proposed, from effective field 
theory descriptions IHZ] and directly QCD-motivated models [HHl to phenomenological 
models 69 . In ref. TO" , a Lagrangian motivated by chiral perturbation theory has been 
fitted to experiment. In paper IV, the tttt transitions from excited QQ states has been 
described by a derivative Qmr coupling, mediated by a heavy scalar resonance. 

5.1 The width for a tttt transition 

The TTTT width of an excited heavy flavor meson is of the form 



where ka and are the four- momenta of the emitted pions, Pi and Pf are those of 
the initial and final state quarkonia, while iOa and LUt denote the energies of the emitted 
pions, respectively. The factors M/E are normalization factors for the heavy meson 
states similar to those employed in ref. |56| . 



r. 
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Evaluation of the above expression leads to the following form for the differential width 
of a TTTT transition, 



1 1 



4 (27r) 



If 



dqq 



(5.2) 



dz 



Mf 



where the variable z is defined by Q • q 
final state quarkonium are given by 



^a{q. z) {Qf + f ) + LUbiq, z) {Qf - f ) Ef{q) '^^'1 ' 

Qqz, and the energies of the pions and the 



= 

LUb = 

Ef = 



q^l^-Qfqz, 



qy^ + Qfqz, 



Mf. 



(5.3) 
(5.4) 
(5.5) 



In the above expressions, the relative momentum of the emitted pions has been fixed by 
the delta functions in eq. 1)5. and can be expressed as 



{Ef^M,)^-{Aml + q^){Ef^M,f 
A{Ef - M,Y - Aq^z^ 



(5.6) 



The integration limit qf corresponds to the maximal momentum of any one of the final 
state particles, e.g. the final state quarkonium. Thus qf can be calculated as the q-value 
of a transition A' — > AX, where X is a particle with mass Mx = 2771-^. 

The above formalism is adapted for computation of a tttt width using a hadronic matrix 
element Tfi. However, since experiments generally measure the invariant mass y^s^, 
then it is useful to define a dimensionless variable x, 



2to^ 



AM 



(5.7) 



where TOtttt denotes the invariant mass ^/s^ of the two-pion system and AM = Mi 
Mf — 2m7r. The relation between q {— \q\) and m^^ may then be obtained as |S3 



I , _ [Mf - (Mf + m^^)^] [Mf - {Mf ~ m^^f] 
" AMf 

From this relation, the transformation Jacobian may be obtained as 



(5.. 



dq 
dx 



AM 



Mf - {Mf+m^^) 



2 1 



AMf 



{Mf-m^^) 



Mf 



{Mf - m^Tr) 



21 



AMf 



(M/+TO^^)| , 

(5.9) 

where |q| is given in terms of m^rTr by eq. (|5.8() . The width for a tttt transition may thus 
be calculated as 



9/ 



dq- 



dq 



dx 



dJ\rR_dq 
dq dx 



(5.10) 



In the above equation, the latter form turns out to be the most useful, since the ex- 
perimental TTTT spectra are presented either in terms of T~^dT/dx or dV /dx versus x. 
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5.2 TTTT transitions in heavy-light mesons 



The emission of two pions from the Ught constituent quarks is modeled in paper III in 
terms of the conventional chiral pion-quark pseudovector coupling model, augmented by 
a pointlike Weinberg- Tomozawa term. The only free parameter in this model is the axial 
coupling of the light constituent quarks, as the quark masses and other parameters 
of the Hamiltonian model were fitted to the D meson spectrum. In addition to the 
single quark amplitude for two-pion emission, the exchange current contribution to the 
Weinberg- Tomozawa interaction was also considered, and found to interfere destructively 
with the single quark amplitude. 

The pion-quark pseudovector coupling H4.1|l gives rise to Born and crossed Born ampli- 
tudes of conventional form for the emission of two pions from an interacting constituent 
quark. The chiral model for the tttt emission amplitude is completed by the Weinberg- 
Tomozawa (WT) interaction, which is described by the Lagrangian 



£wT = ~7~H ^1 '^M '^a TTa X SpTTa 



(5.11) 



The general isospin decomposition of the tttt emission amplitude for constituent quarks 
is, in analogy with that for nucleons. 



T = S„j,r 



\[n,Ta]T , 



(5.12) 



while the general expression for the amplitudes T+ and T is 

r± =u(p') (A±-z7^Q^i?±)u(p). 



(5.13) 



In these expressions, Q denotes the relative four-momentum Q = (kb — ka)/2 of the emit- 
ted pions. In this notation the Born, crossed Born and Weinberg- Tomozawa amplitudes 
are, respectively 



WT 



A 
A 



7t^Qfj. - '^irriq - 4m; 



ItiQ^L + 2imq - Am 



'ba 



1 [ 1 

2 Vb,Ta\ 



(5.14) 

(5.15) 
(5.16) 



Comparison of these amplitudes with eq. H5.13|l yields the desired expressions for the 
sub-amplitudes and , which are of the form 



A+ = 

A- ^ 

B+ = 

B^ = 



q \ 2 

9a 



0, 



5^ 



2/. 

A 
2U 



4m„, 



Ami 



2 + 4to^ 



1 



(5.17) 
(5.18) 
(5.19) 

(5.20) 
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Here the identities — —s and = —u, where s and u are the invariant Mandelstam 
variables, have been used. These results for the A and B amplitudes are formally equiv- 
alent to the corresponding results for the two-pion emission amplitude for nucleons j64| . 
Note that in eq. H5.2()|l . the contribution from the Weinberg- Tomozawa interaction tends 
to cancel the constant term in the B~ amplitude that arises from the Born terms. If the 
axial coupling constant is taken to equal 1, then this cancellation is exact. 

For calculational purposes, eq. H5.13|l was split in paper III into spin-independent and 
spin-dependent parts according to 

^a^ +i(Ty- f3^. (5.21) 

The spin and isospin summed squared amplitude for a given tttt transition is then ex- 
pressed in the form 

|rp = \T\l^ + \T\l^+\T\l^ + \T\l., (5.22) 

from which the tttt width is obtained by insertion into eq. H5.2() . The explicit expressions 
for the above squared amplitudes, given in paper III, are different for each transition 
because of the summation over spin and isospin. The numerical results for the above 
squared amplitudes are presented in Table 

It is instructive for the determination of the two-quark contribution to the Weinberg- 
Tomozawa interaction to write the Lagrangian H5.11(l in the form of a current-current 
coupling, 

-CwT ^TT2 ^ df^T^a, (5.23) 

where V^a = iipq 7^ ipq Ta is the isovector current of the light constituent quark and 
TTa X dp^TTa IS the Current of the tttt system. Given this expression, it becomes natural to 
describe the irreducible two-quark contribution to the tttt production operator by means 
of two-quark interaction current contributions to the isovector current V^q. 

In the non-relativistic limit, the spatial term in the isovector current V^ta — (Va,iVba) 
of the light constituent quark takes the form 



Va 



2mq 



(5.24) 



where niq is the light constituent quark mass and Pq and the initial and final quark 
momenta, respectively. The expressions for the exchange current contributions for the 
scalar confining and OGE interactions have been calculated in ref. fi9) . and may be 
expressed as 



^Co„f ^ -YsMv,, (5.25) 



-q 



yOge _ 



2ml 



Mr 



{p'q +Pq+ if^Q X '^2) + icTq X k2 



(5.26) 



In paper III, the above operators contribute only to \T\'^_ , and were found to reduce the 
widths for vrvr transitions. However, it was also found that the nonrelativistic treatment 
of the two-quark operators is somewhat unrealistic because of the small mass of the light 
constituent quark. An approximate relativistic treatment akin to that employed for the 
single pion transitions in paper II was therefore used in paper III. The results obtained 
upon addition of the two-quark operators H5.25|l and 1)5. 26() are shown in Table 
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Table 5.1: Numerical results from paper III with single-quark amplitudes and relativistic 
matrix elements for the tttt widths of the spin triplet D2 , Di , Dq and the spin singlet 
Z?* mesons. The individual results from the spin independent (a^) and spin dependent 
amplitudes are shown for g\ — \. 



Transition 


ml. 


|y|2 

1 1 OC 


\T 


2 

/3- 




\T 


I. 




= 1 


g\ = 0.87 


Dl 


— > D'tt-k 


0.896 


2.864 


5.06 


10" 


3 


1.17 


10" 


-3 


3.77 MeV 


2.39 MeV 


Dl 


Dtttv 






6.45 


10" 


2 


7.51 


10" 


-3 


0.07 MeV 


0.05 MeV 


Di 


D*TVTV 


0.367 


1.291 


2.33 


10" 


3 


7.85 


10" 


-4 


1.66 MeV 


1.05 MeV 


Dl 


D*Tvn 






6.54 


10" 


4 


3.30 


10" 


-4 


~ MeV 


~ MeV 


Dt 


— > Dtttv 


2.749 


7.915 














10.7 MeV 


6.80 MeV 


Do 


— ^ D*TVK 


0.020 


0.098 














0.12 MeV 


0.07 MeV 


D'o 


Dtttv 






1.43 


10" 


2 


2.86 


10" 


-3 


0.02 MeV 


0.01 MeV 



Table 5.2: Numerical results from paper III for the most important tttt transitions, upon 
addition of the two-quark contributions associated with the scalar confining and OGE 
interactions to the Weinberg- Tomozawa Lagrangian, for = 1. 







\T\l- 




Total 


Transition 


Rel 


+Conf 


+OGE 


ffl = l 


g\ = 0.87 


D2 ^ D'Tin 


2.864 


2.377 


2.144 


3.05 MeV 


1.82 MeV 


Dl ^ D'Tvn 


1.291 


1.076 


0.974 


1.34 MeV 


0.80 MeV 


Dl Dtttv 


7.915 


6.535 


5.872 


8.62 MeV 


5.17 MeV 



The results for the tttt widths of the positive parity charm mesons with L = 1 in Table I^TI 
reveal a strong sensitivity to the value of g\ as well as the large liyperfine splittings in the 
empirical D meson spectrum. Consequently, some transitions are kinematically favored, 
while others, in particular Dq 0*1:1:^ are strongly inhibited by the small phase space 
available. The tttt widths of the L = 1 D mesons are, therefore, also very sensitive to the 
spin-orbit structure of the Qq interaction, which was also found to be the case for the 
analogous single pion transitions in paper II. In the absence of empirical information and 
definite QCD lattice calculations the energies of the as yet undiscovered Dq and Dl 
mesons were in paper III taken to equal those predicted by the calculation of ref . 38 . 

The predicted total widths of the L = 1 D mesons may be obtained by adding the 
calculated tttt widths in Table 15.11 to those for the tt transitions obtained in paper II. 
With (7^ = 1 the total calculated tttt width in the single quark approximation of the 
1)2(2460) is 3.8 MeV and that of the L»i(2420) is 1.7 MeV. Such an addition brings the 
total calculated width of the (2460) to 19.5 MeV, which is well within the uncertainty 
margin of the empirical value 25j:y MeV In the case of the Di(2420) meson, the 

total calculated width for tt and tttt emission comes to 15.3 MeV, which is close to the 
empirical uncertainty margin of the total width 18.9^3 5 MeV [HI]. However, reduction 
of the value for g\ to 0.87 brings the calculated widths for tt and tttt emission somewhat 
below the empirical values for the total widths. Likewise, employment of the two-quark 
contributions associated with the Qq interaction has the effect of reducing the calculated 
TTTT widths, as shown in Table I5T^ 
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5.3 TTTT transitions in heavy quarkonia 

Theoretical work on the tttt decays of excited heavy quarkonia has demonstrated that the 
empirical energy spectra of the emitted tttt system demands that the pions be derivatively 
coupled to the heavy quarkonium states. This is consistent with the role of the pions 
as Goldstone bosons of the spontaneously broken approximate chiral symmetry of QCD. 
Most models |H3 |HS| IMI IZ^ have dealt with the coupling of two-pions to the heavy 
meson as a whole rather than to its constituent quarks. The satisfactory description 
obtained suggests that the decay amplitude Tfi at the quark level should be a smooth 
function of the tttt momentum q, which is dominated by single-quark mechanisms for tttt 
emission. However, the pion rescattering or pion exchange term that appears naturally 
as a consequence of the coupling of two-pions to constituent quarks has been found, 
in paper IV, to be dominant since it is not suppressed by the orthogonality of the 
quarkonium wave functions. 

It was shown in paper IV that an unrealistically large pion rescattering contribution 
may be avoided if the Qtttt vertex involves an intermediate, fairly light and broad a 
meson, in line with the phenomenological resonance model of ref. |fi9j . An intermediate 
a meson suppresses the contributions from pion rescattering mechanisms, while single 
quark amplitudes are but slightly affected. Together with a relativistic treatment of 
the single quark amplitude, this suppression of the pion rescattering contribution has 
been shown in paper IV to reproduce the expected smooth behavior of the transition 
amplitude. 



Q 



kb 





Figure 5.1: Tree-level Feynman diagrams for the emission of two-pions by heavy con- 
stituent quarks. The two upper diagrams correspond to the single-quark amplitudes 
Tq and Tq respectively, while the two lower diagrams describe the pion rescattering 
amplitudes Tox and T^xc which involve a pion exchange between the heavy constituent 
quarks. 
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If the coupling of the pions to the constituent quark does not involve derivatives of the 
pion field, then agreement with experiment is excluded for the pion invariant mass distri- 
butions in the decays T' ^ T tttt and ip' — *■ J/-0 tttt [HHUBni- Derivative couphngs for the 
pions are also consistent with the role of pions as Goldstone bosons of the spontaneously 
broken approximate chiral symmetry of QCD. The effctive Qtttt interaction Lagrangian 
is therefore expected to have the form 

t^QTlTT = 47rA IpQ df^TTa SpTTa ipQ, (5-27) 

where A is a coupling constant of dimension [MeV]"'^ and ipQj'tpQ denote the heavy 
quark spinors. The total tree- level amplitude for tttt emission can then be expressed in 
terms of single quark and pion rescattering terms, as illustrated in Fig. 15.11 The isospin 
dependence of the dipion-quark coupling then implies that 

\T\l^ = 2\T\l^,^+\T\lo,o. (5.28) 

As a consequence the width for emission of charged pions should be twice that for emission 
of neutral pions, which is in fair agreement with what is found experimentally |32l Ififtj . 

The effect of the strong interaction in the tttt system may be approximately accounted 
for by inclusion of an intermediate scalar meson (cr or glueball) resonance in the vertex. 
This is brought about by modification of the coupling constant A with a relativistic scalar 
meson propagator of the Breit-Wigner type: 



Here denotes the pole position mcr — iTcr/2, and q the four-momentum, of the effective 
scalar (cr) meson resonance. The a resonance appears by infinite iteration of the four- 
pion vertex in the isospin spin channel. Therefore, as pointed out in ref. 69 , it is 
natural to describe the strongly interacting tttt state by a broad a pole rather than by the 
driving term (4-pion vertex) alone. When the modification (|5.29|) is taken into account, 
the expression for the single-quark amplitude becomes 



Tiq — — 16 ttA 



\NQ + q^ + Tl/A 



(K + ^bf - q') 



A^iq. (5.30) 



The nonrelativistic approximation for the matrix element Mi^ is unreliable, as the radial 
5'-wave quarkonium wave functions are orthogonal. A relativistic form for the matrix 
element TMiq, where P = |P|, q= \q\ and P ■ q = Pqv, may be obtained as 



M\°^ = - / dr'r'uf{r') / drrui{r) / dP P^ / dva{P,v,q) 
Jo Jo Jo J-1 





where a{P, v, q) is a factor which includes the quark spinors in the coupling (|5.27|) . 
(P \ I {E + Mq){E' + Mq) ^_P^_-qV^^\ 
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In principle the quark spinors in the coupling (|5.27|) also contain a spin dependent part 
that is proportional to q x P. For the present purposes that contribution turns out 
to be very tiny and may be safely neglected. It has been shown in paper IV that the 
relativistic modifications to the single-quark matrix element increase the magnitude of 
the single quark amplitudes for tttt decay. 



Figure 5.2: Modeling of the pion rescat- 
tering term in Fig. 15.11 in terms of inter- 
mediate a mesons. The four-momenta of 
the a mesons are defined as ki = —k — ka 
and ^2 = k — kt- The crossed rescatter- 
ing diagram in Fig. 15.11 can be obtained 
by interchanging ka and fcf,. The physi- 
cally reasonable approximations fci ~ — fc, 
fc2 ~ k, |fcO| « \kl\ « {uja + ^b)/'2- and 
fco ~ {uj})—LDa)/2 were made in paper IV, to 
allow for a simpler treatment of the triple 
propagators in the pion rescattering ampli- 
tudes. 



Q 




Q 



The pion rescattering amplitude that corresponds to Fig. I5.2l mav be expressed as 



+ M2 + r2/4)(fc2 + m2)(fc2 + M2 + r2/4) 



, (5.33) 



where the momenta are defined as in Fig. 15. 21 Upon simplification of the triple propagator 
according to the recipe of Fig. 15.21 the crossed pion rescattering diagram gives an extra 
factor of 2, yielding the expression 



T2q = -128 7r2A2 



1 fl 



+ 



q^z^ 2 



3^4 -Q}) [-^ci - A\M,2 - M,3)] (5.34) 



where the term which is proportional to the matrix element Mei represents an amplitude 
for D-wave tttt emission. The matrix elements in eq. ()5.34ll may, in the non-relativistic 
approximation, be expressed as 



druf{r)ui{r)jo{Qfr) 
druf{r)ui{r)jo{Qfr) 
druf{r)ui{r)jo{Qfr) 
dr Uf{r)ui{r)j2{Qfr) 





r2/4)2 




{Ml + 


r2/4)2 


47r(X2 




{Ml + 


r2/4)2 


47r(X2 


_^2)2 


{Ml + 


r2/4)2 



-Xr 



47r(X2 - A^y 



2X 

XY^{Xr) 
F2{r). 



2X ^ 



(5.35) 
(5.36) 
(5.37) 
(5.38) 



In the above matrix elements, X is defined as X — \J + rl./A— {uja + '-^hY l'^-, while 
yo(r) denotes the Yukawa function e^^ /r. Note that when the value of k^ exceeds m^, the 
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analytic continuation A ^ —i \/kQ ~ (71| is employed for the matrix element 
Further, Uf(r) and Ui(r) denote the reduced radial wave functions for the final and initial 
state heavy quarkonia, respectively. The function F2{r), which is defined in paper IV, 
is closely related to, and in the limit oo actually reduces to, a Yukawa Y2 func- 

tion j71| . It turns out that the matrix element (|5.38|l is numerically quite insignificant, 
because of the strong suppression caused by the j2 function for small values of Qr. Also, 
the smallness of as compared with k precludes any terms proportional to /cq or 
from playing any major role. 

In view of the many approximations involved in the above treatment of the pion rescat- 
tering terms, a check against an unapproximated calculation is desirable. This is possible 
since the triple propagator in eq. (|5.33l) may also be considered without any approxima- 
tion in fci and k2, at the price of numerically much more cumbersome expressions. By 
means of the Feynman parameterization 

^ ^2 f dxx [ dy— , ^ V ^ ,^ ^, (5.39) 



ABC Jo Jo '^[A{l-x) + Bxy + Cx{l~y)]3 
the two-quark amplitudes of Fig. 15.21 mav be cast in the form 

- -{8nX?{lK-Q} 



[ dxx [ dy{Mi-A^Mii] 
Jo J a 

/ dxx I dy ( — — {1 — 2x + xy) — Q^f{l — xy) + qQ fz{l 

Mil - 2x + xy) + Q}{1 - xy) -f qQfZxil - y)^ 



u2 



-ujaUJbko / dxx dy Mu \ +Tcxc, (5.40) 
Jo Jo 



where the matrix elements are given by 



Ml = 2 / druf{r)ui{r) 



/o 




— Ar 



2x + xyY + Qjx'^y'^ + qQfz{l - 2x + xy)xy , (5.41) 



Mil = 21 druf{r)u,{r)^^^^±^^e-^^{rA + l) 



/o 




(5.42) 

Here the term proportional to k^ is again only of minor importance. Note that in order 
to obtain the contribution Texc to eq. (|5.4Q(I . it is necessary to make the substitution 
ka kjj, which implies uja <-> lojj and Qf ^ ~Qf- In the above matrix elements, the 
quantity A is now defined a,s A — ^ vri^ — Kq and involves an effective mass and an 
energy transfer variable Kq- These are abbreviations of the following expressions: 
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ml = (^Ml + ^^{l^xy)^mlx{2{\-x){\-y)-xy^) (5.43) 

+ 2{^-Q}^ x{l-x){l-y), 

Ko = u;ail-x)-uJbxil-y) + ko, (5.44) 

where and ko are defined as TOq- — and {ujh — uja)/'^, respectively. Numerical 

comparison of the above formalism with the approximate model of paper IV indicates 
that eq. H5.34|l is accurate to within ~ 3%. 

All pion rescattering matrix elements have, in paper IV, been considered in the non- 
relativistic limit, even though it was noted that that limit is not realistic in the case of 
the single quark amplitudes. This treatment is expected to be permissible in the case of 
two-quark amplitudes, since the Yukawa functions from the propagators of the exchanged 
pions and a mesons cancel the orthogonality of the radial S'-wave quarkonium wave func- 
tions. Relativistic effects thus constitute only a correction to the pion rescattering matrix 
elements, which is expected to be rather small because of the large constituent masses 
of the charm and bottom quarks. 



5.4 The transitions T' ^ Ttttt and i/j' J/i/jtttt 

If a tr meson lighter than about 1 GeV is employed, together with a relativistic treatment 
of the single quark amplitudes, then the effects of pion rescattering diagrams may be 
reduced, and they may become subdominant as compared to the single quark amplitudes, 
which allows for agreement with experiment. The results of paper IV indicate that a a 
mass of ^ 500 MeV gives a favorable description of the present experimental data on the 
TTTT transitions from the 25* states of heavy quarkonia. 

The calculated widths and tttt energy distributions were obtained by simultaneously 
optimizing the results for T' — )■ Ttt+tt^ and ip' — )■ J/ipTr~^Tr~. The best results from 
paper IV, which yielded the cr meson parameters ttIc — 450 MeV and = 550 MeV, 
are shown in Table I^T^ and Fig. 15.31 A coupling constant close to A = — 0.02 fm'^ was 
found to provide an optimal description of the tttt widths and energy spectra, where the 
sensitivity to the sign of A is due to the consideration of pion rescattering amplitudes. 
The heavy quark masses correspond to those of the QQ Hamiltonian model in paper IV. 



Table 5.3: Experimental data [S] and calculated widths for tttt transitions from the T' 
and V' states, for A = -0.02 fm^, m„ = 450 MeV and T„ = 550 MeV, together with 
experimental widths, branching fractions, and maximal momenta for each transition. 



Transition 


Ttot 


% of Ttot 


Tcxp 


Tcalc 


^max 


T' ^ Ttt+tt- 


44 ± 7 keV 


18.8 ± 0.6 % 


8.3 ± 1.3 keV 


5.89 keV 


475 MeV 


T' ^ T 7r°7r" 




9.0 ± 0.8 % 


4.0 ± 0.8 kcV 


3.07 keV 


480 MeV 


tp' J/lp TT^TV^ 


277 ± 31 keV 


31.0 ± 2.8 % 


86 ± 12 keV 


53.5 keV 


477 MeV 


iP' J/iP ttV 




18.2 ± 2.3 % 


50 ± 10 keV 


27.8 keV 


481 MeV 
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Two-pion decay distribution - Y(2S) -> Y(1 S)7tV 
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Two-pion decay distribution - W(2S) -> S)jr"^Tt' 





— • — Experiment 

Predicted distribution (dimensioniess units) 




2.5 






2 




t/ i - 


1.5 




i 




1 








0.5 



..-4 







0.2 0.4 0.6 0.8 1 



Figure 5.3: Comparison of calculated and experimental |66| tttt energy distributions 
for T' ^ Ttt+tt- and ^p' J/i/'Tr+vr-, for ^ 450 MeV, r<,_= 550 MeV and 
A = —0.02 fm'^. The calculated width for tt+tt^ decay is 5.89 keV for bb and 53.5 keV for 
cc. The scaled tttt invariant mass x is defined in eq. H5.7|l . 

The results presented in Fig. 15. 31 indicate that the shapes of the experimental tttt spectra 
for bb and cc are slightly different. In particular, the peak at high x appears somewhat 
lower for bb, while the tail at low x is more pronounced for bb. It was therefore noted 
in ref. [MI that the resonance model of ref. IMI cannot be simultaneously fitted to both 
the bb and cc data. This is because the shape of the tttt energy distribution is rather 
insensitive to the properties of the a meson when only amplitudes of the single quark 
type are considered. Thus widely different masses and widths of the a meson have to be 
employed to fit the empirical tttt energy spectra in the cc and bb systems. 
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However, if pion rescattering amplitudes are considered as well, then a a mass higher 
than about 450 MeV was shown in paper IV to lead to unrealistically large pion rescat- 
tering contributions, as the tttt energy spectrum then begins to develop a second peak at 
low X. As the pion rescattering effects are of a short-ranged character, then it turns out 
that they are significant only for the tttt spectrum in the bottomonium system. While 
the pion rescattering contribution is seen from Fig. 15.31 to account for the qualitative 
differences between the tttt spectra for hh and cc, it nevertheless appears to be some- 
what overpredicted. This problem can be traced to the nonrelativistic treatment of the 
pion rescattering contribution, and may be alleviated if relativistic matrix elements are 
employed, as discussed in paper IV. 

It may be seen from Table 15.31 that the tttt widths of the 66 system are somewhat un- 
derpredicted. This is because the presence of the pion rescattering terms precludes an 
independent fit of the width and the tttt spectrum. A larger value of A, which would 
improve the tttt widths, would then worsen the description of the tttt spectrum presented 
in Fig. 15.31 However, in the case of the transition V'' J I'^P tt+tt", it may actually be 
desirable to employ a 20 — 30 % larger value of A, as was also suggested in ref. |70| . 



Two-pion decay distribution - Y(3S) -> Y(1S)7t^7t" 
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Figure 5.4: The empirical double-peaked T(3S') Tttt: spectrum, fitted in paper IV 
by the parameters m„ = 1400 MeV, F^ = 100 MeV, A = 2.7 • 10"^, yielding a width of 
F^+^- = 1.07 keV. 



An outstanding problem, lately verified by experimental reanalysis "72' , is the double- 
peaked structure of the empirical T(3S') — > Ttttt spectrum, which cannot be explained 
by models dominated by single-quark amplitudes, such as the one employed in paper IV. 
However, it was also shown that single-quark -I- pion rescattering models do in fact 
have sufficient freedom to accommodate a double-peaked tttt spectrum, as may be seen 
from Fig. 15.41 As a much heavier scalar meson is employed, the contributions from the 
pion rescattering and single quark amplitudes are of equal magnitude. Incidentally, the 
best fit parameters fall within the range of the empirical scalar resonances /o(1370) and 
/o(1500) both of which possess a strong coupling to tttt. 



Chapter 6 

Conclusions 



The electromagnetic and pionic transitions in the heavy QQ and heavy-Hght Qq mesons 
have been calculated within the framework of the covariant BSLT equation, with the 
assumption that the quark-antiquark interaction can be modeled in terms of a long- 
range confining interaction and a short-ranged OGE or instanton induced interaction. 
It has been demonstrated that a reasonable description of the empirical heavy meson 
spectra can be achieved within such an approach, which also yields reasonable values 
for the confining string tension, constituent quark masses and the parameters Aqcd and 
nig which control the strong coupling ois. However, the question of the effective Lorentz 
structure of the quark-antiquark interaction cannot be answered by the quarkonium 
spectra alone, since there are many models which produce equally satisfactory spectra 
using different assumptions for the Lorentz structure of the effective confining interaction. 
Numerical lattice QCD is unrevealing in this the different components of the Qq 

interaction that can be measured on the lattice may also be well fitted by different 
assumptions for the effective confining interaction. 

In this situation, a study of electromagnetic and pionic transitions is instructive, as two- 
quark or negative energy contributions have been shown, within the framework of e.g. 
the Schrodinger and Gross equations, or the instantaneous approximation to the Bethe- 
Salpctcr equation, to be significant for several of these transitions. As the computation 
of these effects requires an explicit assumption of the Lorentz structure of the confining 
interaction, then it is possible that the question may be settled in the future when 
adequate experimental data on electromagnetic and pionic transitions is available. For 
the time being, the Ml transitions in the charmonium system are the most instructive, 
as the two-quark effects for those transitions have been shown to be large. Furthermore, 
as the OGE interaction docs not contribute any two-quark operator for Ml transitions 
in equal-mass quarkonia, then the transition J /ip rjc^ may provide direct insight into 
the Lorentz structure of the confining interaction. Such a calculation has been described 
in this thesis, where it is found that the two-quark operator associated with a scalar 
confining interaction can explain the observed width of about 1 kcV. 

The theoretical importance of the Ml transition J/tp ^ rjc'j suggests that a new and 
more accurate measurement of the width for that transition should be performed as 

soon as possible. The experimental situation is similar to that for the El transitions 
XcJ J/i^l^ which were thought to be overpredicted by most model calculations for a 
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long time, until the issue was resolved by the new data on these El transitions reported in 
the latest edition of the PDG j^Hj. Thus, in order to avoid prolonged speculation, a new 
and independent measurement of J/'0 — > rjcj is desirable. More empirical information is 
also needed on the total widths of the bb states and their El branching fractions, as the 
El widths are very sensitive to the QQ wave functions, if not to the Lorentz structure of 
the quark-antiquark interaction. Progress has recently been made for the El transitions 
in the bb and cc systems as well, since models that employ nonperturbative hyperfine 
interactions have become available. It was found in the calculation reported in this 
thesis that many of the El transitions in the bottomonium system cannot be accurately 
modeled with spin-averaged wave functions. 

A realistic description of the heavy-light Qq mesons, most notably the charmed D meson, 
presents much more serious theoretical challenges for a number of reasons, most notably 
the uncertain composition of the Qq interaction, the relativistic nature of the light con- 
stituent quark and the limited empirical knowledge of the excitation spectra. It has 
nevertheless been demonstrated in this thesis that the Ml transitions provide an instruc- 
tive test for the Lorentz structure of the Qq interaction Hamiltonian. Promising results 
have been obtained for a Hamiltonian with scalar confining and vector OGE components, 
possibly augmented by an instanton induced interaction. The single pion transitions of 
the D meson are likewise instructive in this sense since the two-quark contributions to the 
axial charge component of the transition amplitude are large. Consequently, transitions 
that involve 5'- wave pion emission are predicted to be suppressed, an effect which has 
also been observed within the framework of the Gross equation. 

The recently measured total width of the D* vector meson has been shown to provide 
useful and constraining information on the pion-quark axial coupling constant g^. Also, 
the flavor symmetry violating D* — > Dsir^ transition can not only provide constraining 
information on the magnitude of 7r° — ?7 mixing, but can also be used to estimate the 
magnitude of the 77-quark, and in particular, the yy-nucleon coupling fnNN- In addition to 
single pion emission, tttt transitions have also been found to contribute significantly to the 
predicted total widths of the L = 1 D mesons, a situation which is similar to that in the 
well explored K meson spectrum. These tttt transitions were found to be very sensitive, 
both to the value of g\, as well as to the hyperfine splittings in the D meson spectrum. 
In all likelihood, the tttt transitions will be found to contribute several MeV to the total 
widths for strong decay. The empirically strong tttt transitions in the charmonium and 
bottomonium systems have been investigated by means of a phenomenological model, 
where the Qtttt interaction is mediated by a broad scalar a meson. It has been shown 
that such a model can explain several features of the empirically observed tttt transitions 
in heavy quarkonia, although a completely satisfactory description was not achieved. 

An important conclusion reached in this thesis is that the most instructive test for a 
given QQ or Qq interaction Hamiltonian is not the excitation spectrum but rather the 
radiative Ml transitions between the ground state vector and pseudoscalar mesons. Once 
the question of hyperfine splittings and total widths has been settled by experiment, then 
the pionic transitions from the L = 1 D mesons will provide a similar testing ground. 
So far, a pure scalar confining interaction has passed the above tests, although other 
conceivable forms have not been systematically ruled out. 



Svensksprakigt sammandrag 



Denna avhandling presenterar en utrakning av elektromagnetiska och starka overgangar 
i mesoner med en (Qq) cller tva (QQ) tiuiga kvarkar. Dessa niesoncr har beskriv- 
its med hjalp av den kovarianta Blankenbecler-Sugar (BSLT) ekvationen, under an- 
tagandet att vaxelverkningen mellan kvarkarna kan beskrivas som summan av en ef- 
fektiv fjattrandc viixclvcrkning mod lang rackvidd och cn kort-rackviddskomponent, 
vilken ger upphov till hyperfinstruktur i mesonernas excitationsspektrum. En dylik 
vxelverkning ger en acceptabel beskrivning av det empiriska mass-spektret, ifall kort- 
rackviddskomponcntcn bcskrivs mod hjlp av storningstcorctiskt glnon-utbytc (OGE) cUcr 
en instanton-inducerad vaxelverkning. Emellertid kan mass-spektret inte ensamt avgora 
den relativistiska Lorentz-strukturen hos den fjattrande vaxelverkningen, eftersom flera 
olika ansatser ger ekvivalenta beskrivningar av det empiriska spektret. 

I denna situation ar de elektromagnetiska och starka overgangarna av stor betydclsc, 
eftersom dessa har visats vara kansliga for negativ-energi bidrag till overgangsamplituden, 
vilka i sin tur beror uttryckligcn pa vaxelverkningens Lorentz-struktur. Dessa effekter 
har, till dags date, pavisats inom ctt antal tcorctiska beskrivningar, diiribland Sclirodinger- 
och Gross-ekvationerna. Eftersom utriikningen av dessa bidrag till sonderfallsvidderna 
kraver en ansats for kvark-antikvark vaxelverkningens Lorentz-struktur, kan de starka 
och elektromagnetiska sondcrfallcn utgora ctt test for olika modeller fr den fjattrande 
vaxelverkningen. For ogonblicket ar den magnetiska (Ml) dipolovergangen J/ijj —^ 
rjcj av storsta betydelse, eftersom den experimentellt uppatta vidden pa ~ 1 keV 
bcskrivs daligt av den ickc-rclativistiska kvarkmodcUcn, vilken Icdcr till en tredubbel 
overuppskattning av detta resultat. Ett av nyckelresultaten i denna avhandling ar, att 
en skalart kopplad effektiv fjattrande vaxelverkning kan forklara den uppmatta vidden 
pa ^ 1 kcV. Emellertid bor cn ny cxpcrimcntell matning utforas innan en definitiv slut- 
sats kan dras av det ovanstaende resultatet. 

En annan slntsats, presenterad i denna avhandling, angar de elektriska (El) dipolsonderfallen 
i charmonium (cc) och bottomonium (bb). I detta fall visar det sig att Lorentz-struktren 
hos kvark-antikvark vaxelverkningen i de fiesta fall endast har en forsvinnande liten effekt 
pa overgangsamplituden. Daremot ar El overgangarna kansliga for sma effekter i meson- 
ernas vagfunktioner. Darmed kan en realistisk beskrivning av ett flertal sonderfall endast 
uppnas, ifall kvark-antikvark vaxelverkningens hyperfinstruktur behandlas fuUstandigt, 
vilket ar mqjligt ifall mesonerna beskrivs med hjalp av BSLT-ekvationen. 

Jamfort med QQ mesonerna staller en realistisk beskrivning av Qq systemet myckct 
stora krav pa de teoretiska modellerna, eftersom vaxelverkningens form mellan tunga 
och latta kvarker ar osaker, och den latta kvarken i hog grad relativistisk. Dessutom 
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kompliceras situationen ytterligare av den knapphandiga empiriska kunskapen om Qq 

mesonernas mass-spcktrum. Into dcsto mindre har det visats i dcnna avhandling, att 
Ml overgangarna i D mesonerna utgor ett viktigt test for Qq vaxelverkningens Lorentz- 
struktur. Lovande resultat har erhallits for en skalar fjattrande + OGE vaxelverkning, 
mojligen med tillsats av en instanton-inducerad komponent. De starka pion-sondcrfallcn 
inom D mesonerna ar ocksa av stort intresse eftersom negativ-energi bidragen till den 
axiala laddningsamplituden ar stora. En forutsagelse presenterad i denna avhandling 
ar att pion-overgangar drivna av den axiala laddningsamplituden ar starkt forhindrade, 
vilkct nyligcn ocksa har obscrvcrats mod hjalp av Gross-ckvationcn. 

Den nyligen uppmatta totala vidden for den exciterade D* mesonen ar av stor betydelse, 
eftersom den kan fixera vardet pa den axiala kopplingskonstanten g\ for latta kvarkar. 
Dcssutom kan dot smak-symmctri brytandc sondcrfallct D* DsTt'^ gc information om 
styrkan hos rj mesonens koppling till kvarkar och baryoner, forutsatt av storleken hos 
r] — TT^ blandningsvinkeln ar kand. Utover sonderfall, i vilka endast en pion emitteras, 
kan jivcn tva-pion (tttt) sonderfall vara av betydelse i D mesonerna. Dcnna slutsats 
overensstammer med den experimentella situationen i de sara K mesonerna, dar tttt 
sonderfallen ar val uppmatta. I denna avhandling befanns tttt sonderfallen vara mycket 
kansliga bade for vardet pa och det tillgangliga fasrummct. Dot ar darigenom sanno- 
likt, att TTTT sonderfallen utgor flera MeV av de totala vidderna for starkt sonderfall i D 
mesonerna. 

De empiriskt betydande tttt overgangarna i charmonium (cc) och bottomonium (bb) har 
i denna avhandling undcrsokts mod hjalp av en fcnomcnologisk modcll, vilkcn bcskriver 
kvark-pion vaxelverkningen med hjalp av en skalar a resonans. En dylik modell befanns 
ge en god beskrivning av flera egenskaper hos tttt overgangarna, aven om en fullstandigt 
tillfredstallande beskrivning inte uppnaddes. 

Den viktigastc slutsatscn i dcnna avhandling iir att det most bctydclscfulla tcstct for en 
given kvark-antikvark vaxelverkningsmodell ar inte mass-spektret, utan snarare de elek- 
tromagnetiska Ml overgangarna dar mesonernas grundtillstand byter spinn. I framtiden, 
nar hypcrfin-nivacrna och de totala vidderna hos dc cxcitc;radc; D mesonerna iir cixper- 
imentellt kartlagda, kommer dessa att utgora ett ytterligare test for de ovannamnda 
vaxelverkningsmodellerna. Hittills har en skalar fjattrande vaxelverkning klarat dessa 
test, aven om andra tankbara former inte har uteslutits systematiskt. 
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Suomenkielinen tiivistelma 



Tama vaitoskirja kasittelee sahkomagneettisia ja vahvoja siirtymia mesoneissa, jotka 
koostuvat joko kahdcsta raskaasta {QQ) tai yhdcsta raskaasta ja yhdcsta kcvycsta [Qq) 
kvarkista. Naita mesoneja on kuvattu kovariantin Blankenbecler-Sugar (BSLT) yhtalon 
ratkaisuina, olettaen etta kvarkkien valinen vuorovaikutus voidaan esittaa efektiivisen, 
pitkan kantaman kvarkkcja kahlitscvan vuorovaikutukscn ja lyhycn kantanian ylihicno- 
vuorovaikutuksen summana. Mikali lyhyen kantaman vuorovaikutus oletetaan joko 
hairioteoreettiseksi gluonivaihdoksi taikka instantoni-indusoiduksi, on tuloksena laadul- 
taan hyvaksyttava malli mcsonicn kokccUisillc viritysspcktrcillc. Valitcttavasti viritys- 
spektrit eivat yksinaan riita antamaan ratkaisevaa tietoa kahlitsevan vuorovaikutuksen 
suhteellisuusteoreettisesta Lorentz-rakenteesta, silla useampien eri oletusten on todettu 
johtavan samanlaatuisiin spektreihin. 

Sahkomagnccttistcn ja vahvojcn vuorovaikutustcn mcrkitys on tassa tilantccssa crittain 
suuri, silla on osoitettu naiden riippuvan herkasti siirtymaamplitudin negatiivisen 
energian komponenteista. Nama puolestaan riippuvat eksplisiittisesti kvarkkien valisen 
vuorovaikiitiiks(ni Lorcntz-rakcntccsta. Tahan johtopaatokscc!n on paiidytty aikaiscmmin 
m.m. Schrodingerin ja Grossin yhtaloiden kautta. Koska siirtymaamplitudien laskeminen 
vaatii oletuksen kvarkkien valisen vuorovaikutuksen Lorentz-rakenteesta, voivat nama 
siis toiniia tc!stin;i, milla voidaan vcrtailla cri mallicn todenmukaisuutta myos silloin, kun 
viritysspektrien ennustukset ovat degeneroituneet. Talla hetkella magneettisen (Ml) 
dipolisiirtyman J/ip — > t?c7 merkitys on hyvin suuri, silla sen kokeellinen viivaleveys 
^ 1 kcV ci ole sopusoinnussa eparclativistiscn cnnustuksen kanssa, joka on kolmc kertaa 
suurempi. Eras taman teoksen tarkeimmista tuloksista on, etta skalaarinen kahlitseva 
vuorovaikutus tarjoaa mahdoUisen selityksen yllamainituUe kokeelliselle viivaleveydelle. 
Lopullinen johtopaatos vaatii kuitenkin kokeellisen tuloksen varmistamisen uudella mit- 
tauksella. 

Tassa vaitoskirjassa on myos kasitelty raskaiden mesonien sahkoisia (El) dipolisiirtymia, 
jolloin on todettu etta kvarkkien valisen vuorovaikutuksen Lorentz-rakenteella on 
useimmiten vain haviavan pieni vaikutus teoreettisiin viivaleveyksiin. Toisaalta El siir- 
tymat ovat herkkia pienille muutoksille mesonien aaltofunktioissa. Nain oUen onkin tassa 
teoksessa osoitettu, etta useita El siirtymia voidaan kuvata onnistuneesti vain, mikali 
kvarkkien valinen ylihieno-vuorovaikutus otetaan huomioon kayttamatta ensimmaisen 
kertaluvun hairioteoriaa. 

Mikali mesoniin sisaltyy kevyt kvarkki, on realistisen mallin rakentaminen heti huo- 
mattavasti vaikeampaa, koska raskaiden ja kevyiden kvarkkien valisen vuorovaikutuk- 
sen muoto on erittain kyseenalainen. Lisaksi kevyt kvarkki on hyvin relativistinen, ja 
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viritysspektrin kokeellinen tuntemus vahainen. Tasta huolimatta on tassa teoksessa os- 
oitcttu, ctta D mcsonien Ml siirtymat tarjoavat mahdoUisuudcn tutkia Qg-viiorovaikutuksen 
Lorentz-rakennetta, joUoin lupaavia tuloksia on saatu edellamainituilla vuorovaikutus- 
malleilla. D mesonien vahvat pionisiirtymat ovat myos tassa mielessa kiinnostavia, silla 
pionin ja kcvycn kvarkin valisen vuorovaikutukscn aksiaalincn varauskomponentti riip- 
puu vahvasti kvarkkien valisen vuorovaikutuksen Lorentz-rakenteesta. Tassa vaitoskirjassa 
on havaittu, etta aksiaalisesta varauskomponentista riippuvat siirtymat lienevat voimakkaasti 
cstyneita. Samankaltaiseen johtopaatokseen oUaan paasty hiljattain myos Grossin yhtalon 
kautta. 

Eksitoituneen D* mesonin viivaleveys, josta hiljattain saatiin ensimmainen kokeellinen 
mittaus, on teoreettisesti hyvin tarkea suure, silla sen avulla voidaan maarata kevyi- 
den kvarkkien aksiaalincn kytkinvakio g^. Lisaksi voidaan makusymmctriaa rikkovan 
D* — > DsTT^ siirtyman avulla tutkia r] mesonin ja kvarkkien (tai baryonien) valista 
vuorovaikutusta. Tama edellyttaa, etta neutraalin pionin ja 77 mesonin valinen sekoi- 
tuskulma on timncttu. Lisaksi on tassa vaitoskirjassa tutkittn kalidcn pionin (tttt) siir- 
tymia, jotka hyvin todennakoisesti ovat merkityksellisia D mesoneissa, mika johtopaatos 
on yhtapitava oudosta K mesonista saadun kokeellisen tiedon kanssa. Tassa teok- 
sessa on osoitcttu D mcsonien tttt siirtymicn olcvan herkkia g\'.Ti numccriscUc arvoUc 
seka kaytettavissa olevalle faasiavaruudelle. Taten tttt siirtymien viivaleveydet ovat to- 
dennakoisesti muutaman MeVm suuruisia. 

Charmonium (cc) ja bottomonium (66) mesonien kokeellisesti merkittavia tttt siirtymia 
on tassa vaitoskirjassa tutkittu fcnomcnologiscUa mallilla, jossa raskaidcn kvarkkien ja 
pionien valista vuorovaikutusta kuvataan skalaarin a resonanssin avulla. Talla mallilla 
saavutettiin tyydyttava, joskaan ei taydellinen, kuvaus raskaiden mesonien tttt siirtymista. 

Taman vaitoskirjan tarkein johtopaatos on, etta sahkomagneettiset Ml siirtymat, joissa 
kvarkkien spin rnuuttuu, rnuodostavat viritysspcktrista riippnmattoman tcstin kvarkkeja 
kahlitsevan vuorovaikutuksen Lorentz-rakenteelle. Tulevaisuudessa, ylihienorakenteen 
ja viivaleveyksien ollessa tunnettuja, tulevat D mesonien Ml ja pionisiirtymat tarjoa- 
rnaan vastaavan tcstin. Talian mcnncssa on skalaarincn yritc kvarkkeja kahlitscvalle 
vuorovaikutukselle osoittautunut menestyksekkaaksi, vaikka muita mahdollisia muotoja 
ei ole jarjestelmallisesti eliminoitu. 
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